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Abstract 



We study a class of homeomorphisms of surfaces collectively known as linked- twist maps. 
We introduce an abstract definition which enables us to give a precise characterisation of 
a property observed by other authors, namely that such maps fall into one of two classes 
termed co- and counter-twisting. We single out three specific linked-twist maps, one 
each on the two-torus, in the plane and on the two-sphere and for each prove a theorem 
concerning its ergodic properties with respect to the invariant Lebesgue measure. 

For the map on the torus we prove that there is an invariant, zero- measure Cantor set 
on which the dynamics arc topologically conjugate to a full shift on the space of symbol 
sequences. Such features are commonly known as topological horseshoes. For the map 
in the plane we prove that there is a set of full measure on which the dynamics are 
measure-theoretically isomorphic to a full shift on the space of symbol sequences. This 
is commonly known as the Bernoulli property and verifies, under certain conditions, a 
conjecture of Wojtkowski's. We introduce the map on the sphere and prove that it too 
has the Bernoulli property. 

We conclude with some conjectures, drawn from our experience, concerning how one 
might extend the results we have for specific linked-twist to the abstract linked-twist 
maps we have defined. 
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1 Introduction 



The work in this thesis can be categorised as dynamical systems or non-linear dynamics. 
This huge field, in broad terms, studies the trajectories of the points which constitute 
some space, given some rule which governs the evolution of that space as time progresses. 
It has strong connections to many of the major fields in pure and applied mathematics, 
to the natural sciences and to engineering. The present work is primarily of a pure- 
mathematical nature and relies heavily upon the results and techniques of ergodic theory. 

Ergodic theory studies dynamical systems with an invariant measure. We discuss 
ergodic theory in greater detail in Section 12.11 Ergodic theory is built upon measure 
theory, itself one of the cornerstones of mathematical analysis. Its influence is felt in 
two crucial ways: it allows us to describe and to prove certain limiting behaviour, which 
provides us with information about the evolution of our dynamical system; and it allows 
us to disregard certain points which evolve in a manner that is atypical and inconvenient 
for us. 

Similarly important is hyperbolicity, which we discuss in Section 12.21 Hyperbolic 
behaviour in our dynamical systems is of critical importance insofar as all of our tech- 
niques for demonstrating ergodic properties rely upon it. In essence (and of course, 
we give rigorous definitions later) hyperbolicity concerns the behaviour of those points 
'close to' some reference point whose evolution we are following. Depending upon the 
direction of the displacement, these nearby points either approach or move away from 
our reference point as we evolve the system, but crucially they do not stay at a fixed 
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distance. This behaviour can lead to initial conditions being perpetually thrown apart 
and back together and result in a mixing of the ambient space. 

In the remainder of this introduction we will introduce the maps that we shall study 
and state the three main theorems we shall prove. We do not do so by the most direct 
route however, preferring first to motivate the concept of hyperbolicity in a simple 
example. This occupies Section 11.11 In Section 11.21 we introduce the reader to the 
linked-twist maps with whose properties this work is concerned. We do this first in an 
abstract setting which enables us highlight what unites them all and classify them in 
an important way. Finally Section [1.31 is divided into three parts, in each of which we 
define a linked-twist map and state a theorem we shall prove for that map. 

Following on from this, the remainder of our thesis is organised as follows. In Chap- 
ter [2] we provide a literature review which is divided into four sections. In Section 12.11 
we discuss ergodic theory, providing the definitions we will need throughout this work, 
in particular of the Bernoulli property. In Section 12.21 we discuss hyperbolicity and de- 
scribe some important results we will use. In Section [2. 31 we survey those results already 
known for the maps we shall study. Lastly in Section 12.41 we shall discuss a number of 
applications which can be modelled by linked-twist maps. 

Chapters [3l m and [5] are where we prove the new results. In each case we define the 
map and state the theorem later in this introduction, then give a detailed breakdown 
of the method at the start of the chapter. In Chapter [3] we show that a linked-twist 
map defined on a subset of has an invariant, zero-measure Cantor set on which the 
dynamics are topologically conjugate to a full shift on N symbols. For further details 
see Section 11.3.11 In Chapter H] we show that a linked-twist map defined on a subset 

jesgu e measure. This verifies 



of the plane has the Bernoulli property on a set of full Lei 

(under certain conditions) a conjecture of Wojtkowski's (|l98Cll ). a precise statement of 
which is postponed until that chapter, where we establish the required notation. We 
give more details in Section ri.3.2[ Finally in Chapter [5] we prove the Bernoulli property 
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for a linked- twist map defined on a subset of S^. We introduce this map in Section ri.3.31 
We conclude in Chapter[6]by analysing the results we have established and discussing 
the strengths and weaknesses of our methods. There are some obvious generalisations 
which suggest themselves as well as some different directions one could take whilst still 
building upon the work we have done, so we consider both. Based on what we have 
learned we feel confident in making some conjectures and we include these here. 

1.1 Motivation 

We begin by describing a system which illustrates hyperbolicity in perhaps the simplest 
non-trivial setting. We will use some of the language of ergodic theory and hyperbolic 
theory to be introduced in Sections 12.11 and 12.21 The reader who is unfamiliar with 
these terms is encouraged to skip forward to these definitions as necessary, although we 
have tried to keep the exposition as elementary as is possible. 

1.1.1 A hyperbolic toral automorphism 

Hyperbolic toral automorphisms are canonical examples of dynamical systems displaying 
hyperbolic behaviour. We describe one here, commonly known as the cat map. More de- 



Katok and Hasselblatt (19951 ) 



tai ls can be found in most dynamical systems text; we recommend[ 

or lBrin and Stuck (2002 ). Given the two-torus = M^\Z^, the cat map is the linear 







diffeomorphism H :T'^ —>-T'^ given by | 

H{x, y) = {x + y,x + 2y) mod 7?. 

We naturally think of as the unit square in the plane with opposing sides identified. 
In Figure 11.11 we illustrate H by first viewing it as a linear map of the plane and then 



^It is perhaps more common to define the map as (x, y) {2x ~\-y,x-\-y) but this is merely a matter 
of personal taste and the results we will quote hold for any hyperbolic toral automorphism. When we 
introduce linked- twist maps on we will wish to emphasise the cat map as a special case, and for this 
purpose our definition is more convenient. 
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seeing how the 'pieces' fit back together on T^. 




Figure 1.1: The 'cat map'. Part (a) illustrates which we represent as the unit square 
in the plane. The shading will help us to illustrate the map. Part (b) shows the image 
of the unit square under H if we consider as a linear map of the plane (i.e. without 
taking the image modulo 7?). Part (c) shows how this image looks upon projection to 

Let us describe, without giving the general definition, what we mean when we say 
that the cat map is hyperbolic. The Jacobian matrix is given by 




and is independent of z G T^. It has distinct real eigenvalues < A_ < 1 < A+ = 1/A_ 
and corresponding eigenvectors v± = (1, A± — 1). Using only elementary linear algebra 
we can draw some simple conclusions about the dynamics of H. 

Suppose that z G and consider the line through z having gradient v-\ we call 
this line the stable manifold of z. It is easily checked that the gradient is irrational 
and so the line extends indefinitely and never self-intersects. Let z' = z + kv^, where 



4 



E M, be on this line. lf| Then H{z') = H{z + kv_) = H{z) + H{kv^) = H{z) + kX_v_, 
i.e. H{z') is in the unstable manifold of H{z). Moreover the distance between the 
points (as measured along the unstable manifold) is smaller by a factor of A_ than the 
corresponding distance between z and z' . 

We can repeat this construction using in place of V- to obtain the unstable 
manifold of z. In this case the distance between points is increased by a factor of A+. 
These facts together show that the cat map is hyperbolic; in fact we can say more than 
this. The picture to have in mind is of the two distinct (in fact, orthogonal) directions 
experiencing stretching and contraction respectively. The constructions we have given 
hold for any z G and the growth rates established hold uniformly at each point, so 
in fact we say H is uniformly hyperbolic or even an Anosov diffeomorphism. 

It transpires that from these few facts one can establish a great deal about the 
dynamics of the cat map. In particular it is ergodic, mixi ng and has the B ernoulli 



Katok et al. (1986 ) which 



property. In Section 12.21 we will describe a theorem due to 
gives sufficient criteria for a map to have all of these properties. One could certainly 
use this theorem to establish them for the cat map; however one would, metaphorically 
speaking, be using a sledgehammer to cr ack a walnut. For mor e elegant ways to prove 
such results we recommend the book of iBrin and Stuck (20021 ). in which many more 



than these three properties are established for hyperbolic toral automorphisms. 



1.2 Abstract linked- twist map theory 

In this section we describe what we will call an abstract linked-twist map. The results 
presented in this thesis are all for specific linked-twist maps and the reader who is 
eager to understand the maps we have studied and the results we have proven can 
safely overlook this section on first reading. We would encourage her to return to this 
material later though, for two reasons. 

•^For the benefit of a cieaner exposition we are not appending 'moduio Z^' to our points. 
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First, this section is our attempt to formalise what precisely it is that the different 
maps on different surfaces that are all referred to as linked-twist maps have in common. 
This is perhaps a simple exercise but nevertheless it serves to draw together the results 
we present. 

Second, perhaps more interestingly, we define a property of linked-twist maps which 
divides them into two classes, namely the co-twisting and counter-twisting classes. The 
distinction can have great implications for the dynamics of otherwise similar maps. 
Other authors have noticed this distinction but have treated it as something which 
must be determined for a given map; conversely we define it for an abstract linked-twist 
map and later prove that a given linked-twist map is either co- or counter- twisting. We 
are grateful to Prof. Robert MacKay for his helpful suggestion, from which this idea 
was born. 

1.2.1 Review section: smooth embeddings 

We begin with a number of definitions from the field of differential geometry. The 
terminology will be necessary in order to define an abstract linked-twist map; the reader 
who is already comfortable with the definition of a smooth manifold and an orientation- 
preservi nq embedding can safely skip these. Our definitions a.re taken from , the excellent 



book of iDo Carmo (1976l). We also recommen d the book of 



review section given by 



Brin and Stuck ^2002 1 



Spivak (197(^ ) or the short 



Definition (Smooth manifold of dimension 2). A smooth manifold is a set S together 
with a family of one-to-one maps (pa '■ Ua ^ S of open sets Ua C into S such that 

1- Uq (kaiUa) = S, and 
2. for each pair a, (3 with 

we have that 
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(a) 4>a^{W) and (l)p^{W) are open sets in M^, and 

(b) (f)'^^ o (f)^ and o are differentiable maps. 

The pair (Ua, 4'a) with p G (paiUa) is called a coordinate system of S around p. The 
image (f>a{Ua) is called a coordinate neighbourhood and if g = 4>a{ua,Va) € S, we say 
that {ua,Va) are the coordinates of q in this coordinate system. 

Definition (Orientable; oriented). A smooth manifold S is called orientable if it is 
possible to cover it with a family of coordinate neighbourhoods in such a way that if 
p G S belongs to two such neighbourhoods then the change of coordinates has positive 
Jacobian. The choice of such a family is called an orientation of S and S is called 
oriented. 

Familiar examples of orientable surfaces include the two-torus and the two-sphere 
S^. Conversely the Mobius strip is not orientable. 

We now extend the notion of a differentiable map in the context of smooth manifolds 
of dimension 2. 

Definition (Differentiable map). Let Si and S2 be smooth manifolds of dimension 2. A 
map f : Si ^ S-2 is differentiable at p E Si if given a parametrization : V G M? ^ S2 
around f{p) there exists a parametrization (f> : U C — -Si around p such that 

f((p{U)) C ipiV) and the map 

V'"^ o / o (/>:[/ c ^ 

is differentiable at (/)'~^{p). The map f is differentiable on Si if it is differentiable at 
every p E Si. 

Definition (Immersion). A differentiable map f : S ^ M.^ of S, a smooth manifold of 
dimension 2, is an immersion if the differential 
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Df,:Tp{S)^Tj^j,){R^) 

is injective for each p G S. 

We can now state the definition of an embedding. 

Definition (Embedding). Let S be a smooth manifold of dimension 2. A differentiahle 
map f : S is an embedding if it is an immersion and a homeomorphism onto its 

image. 

Finally, an embedding is called orientation-preserving if its Jacobian has positive 
determinant, and orientation-reversing otherwise. We illustrate the situation in Fig- 
ure O 



^2 h h 




Figure 1.2: Given a differentiahle map E : C ^ M, the figure shows two possibilities 
for the image of the standard basis of M?, denoted (ei, 62) and shown in part (a), under 
the differential DE^. Parts (b) and (c) show bases (/i,/2) of M^, where fj = DEz{ej) 
for j = 1,2. In part (b) DE^ has preserved the orientation, or 'handedness', of the 
standard basis, as shown by the arrow. The corresponding map : C — > M is called 
orientation-preserving. Conversely in part (c) DE^ reverses the orientation of the basis. 
In this case E : C ^ M is called orientation-reversing. 

1.2.2 Abstract linked-twist maps 

Let be the circle. Without loss of generality we assume a coordinate x G [0, 1] on 
where and 1 are identified. In some situations it will be convenient to use some 
other interval in place of [0, 1]; in that case obvious amendments should be made to our 
definitions. 

Let / = [«0)^i] C M be a closed interval. Moreover we will want I C [0, 1] (or in 
the closed interval we use in place of [0, 1] as the case may be). The Cartesian product 



8 



C = S"^ X / is called a cylinder or an annulus and consists of pairs (x, y) such that 
X G §^ and y £ I. C is an oriented smooth manifold (with boundary) and we identify 
the tangent space T^C at a point z £ C with M?. We give T^C the standard basis 
(61,62), where 61 = (1,0) and 62 = (0, 1) in the usual Cartesian coordinates. 
We define a class of homeomorphisms of C = x /: 

Definition (Twist map; twist function). A twist map T : C ^ C is a map of the form 



T{x,y) = {x + t{y),y), 



where t : I ^ , called a twist function, satisfies the following conditions: 

1. t is continuous on [io^ii] and differentiable on {io,ii), 

2. t{io) = and t{ii) = 1 (or an equivalent condition using a different interval for 

3. dt/dy > on {io,ii)- 

We comment that other authors call T an inteqrable twist map. T preserves area 



(Lebesgue measure) and orientation; see 



Katok and Hasselblatt (19951 ). Two possibil- 



ities for the twist function t are shown in Figure II. 3[ Part (a) of the figure illustrates 
a linear twist (we should properly call this an affine twist, of course), of the kind our 
twist maps will be constructed from. It is defined by 



{y - io) / {k - io) ifyG[io,n], 
otherwise. 



:i.2.i) 



The function is not differentiable at y G {ioi^i}- 

Part (b) shows a smooth (i.e. everywhere differentiable) twist of the kind studied 



by 



Burton and Easton (1980 ). It is defined by a cubic equation in y. Smooth twists 
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Sturman et al. (20061 ) for further 



require a different kind of analysis to that which we s hall conduct and we do not intend 
to discuss them in this thesis; see the original paper or 
details. 

11 / — ll 



t{y) 




0^ 

ti 1 ^^-^ io y ^l 1 

Figure 1.3: Linear and smooth twist functions respectively (recall that and 1 are 
identified in S^). Our twist maps (and hence our linked- twist maps) will be constructed 
from the former. Each introduces a different problem into the analysis; the former 
because of the non-differentiable points and the latter because there is no lower bound 
on the derivative. 



We now define an abstract linked-twist map on a subset i? of a two-dimensional 
smooth manifold M. We do so with reference to the cylinders we will embed in M to 
create R. Later on, when we define the linked-twist maps to be studied in this thesis, 
we more commonly do so directly on i? C M. We introduce an important definition. 

Definition (Transversal embedded cylinders). Consider two embedded cylinders in 
some two-dimensional manifold M, i.e. we have cylinders Ci and diffeomorphisms 
Ei : Ci ^ M for i = 1,2. Suppose that Ei{Ci)nE2{C2) ^ and let Zi G Ci be such that 
Ei{zi) = E2{z2) G M . We will say that such embedded cylinders are transversal if and 
only if the vectors {DEi)z^{ei) and (Z) £"2)22(61); which lie in T^^f^^i)^ = ^^2(22)-^' '^^^ 
themselves transversal in the usual sense (i.e. they form a basis for the tangent space). 

We call the connected region(s) Ei{Ci) R E2{C2) the intersection region(s). For 
examples of pairs of transversal embedded cylinders the reader is encouraged to look 
ahead to Figures 11.41 11.51 and 11.81 We can now define a linked-twist map. 

Definition (Linked-twist map). Let M be a two-dimensional oriented smooth manifold 
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and let Ei : Ci ^ M, for i = 1,2, be a pair of transversal embeddings of cylinders 
Ci = S^ X li into M. Denote R = Ei{C) U E2{C) C M. Let Ti : d ^ d for i = 1,2 
be two twist maps given by Ti{x, y) = (x + ti{y),y) where the twist functions : / — > §^ 
satisfy the conditions in the definition above. 
For i = 1,2 and p £ M define Hi : R ^ R by 



E,oTioE-\p) ifpeE,{C), 
id otherwise, 



where id denotes the identity map. A linked-twist map H : R ^ R is given by the 
composition H = H2 o H{ where j and k are positive integers. 

All linked-twist maps of this form can be categorised as either co- or counter-twisting. 
The definition is as follows: 

Definition (Co-twisting; counter- twisting). Let H be a linked-twist map as above and 
let Ei,E2 : C — > M be the transversal embeddings with which it is defined. If both Ei 
and E2 are orientation-preserving, or both Ei and E2 are orientation-reversing, then 
we say that H is counter-twisting. Conversely if one of Ei,E2 is orientation-preserving 
and the other orientation-reversing, then we say that H is co-twisting. 

We have some comments to make regarding the definition. 

First and foremost it might seem to the reader counter-intuitive to give the definition 
as we have, with the co-twisting systems defined as those where the embedded cylinders 
have different orientations; in fact, in light of our definition, we agree. However there is 
a considerable literature for linked-twist maps and we would like our definition to agree 
with it. The terminology seems to have been introduced by 



Sturman et al. (20061 ) and 



their reasoning can be best understood once we have defined linked-twist maps on the 
torus; we do this in the next section. 

Second, the counter-twisting maps (at least, all of the explicit examples of which we 
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are aware) are more difficult to analyse than the corresponding co-twisting maps. In 
this thesis we deal exclusively with co-twisting linked-twist maps so we do not intend 
to say too much about why this is so, but when we survey the literature in Section 12.31 
we will see that, where corresponding co- and counter- twisting maps can be shown to 

have strong ergodic properties, the criteria are more restrictive in the latter case. 

Third, we will dispense entirely with the other notion introduced bv lSturman et al. (20061 ) 



of CO- and counter-rotating linked-twist maps. This notation was intended to explain 
the relative sense of rotation of the two twist maps acting on the embedded cylinders, 
but leads to the somewhat uncomfortable situation whereby planar linked-twist maps 
(introduced in the next section) are simultaneously co-twisting and counter-rotating or 



vice versa. 



1.3 Definitions and statements of theorems 

We now introduce the linked-twist maps to be studied and state the theorems we shall 
prove. All of these maps fit the abstract definition we have given above, though we shall 
not prove so in every case. In the first two cases it is quite obvious. In the case of the 
third map the embedding uses functions with which the reader may not be familiar so 
we will provide all the details. We shall prove that each map is co- twisting. 

1.3.1 Linked-twist maps on the two-torus 

The simplest linked-twist maps to define and analyse are those on the torus. In this 
section we will define a toral linked-twist map and state a theorem to be proven in 
Chapter [3l We give an overview of results in the literature for toral linked-twist maps 
in Section 12.3.11 In Section 12.4.11 we discuss a situation where toral linked-twist maps 
can be used to model the behaviour of certain physical phenomena. As mentioned above 
we will define the map directly on the torus. 

Let §^ denote the closed unit interval [0, 1] with opposite ends identified. We identify 
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the two torus, denoted with the Cartesian product S-*^ x S^. This gives us two angular 
coordinates {x,y). 

Fix four constants < xq < xi < 1 and < yo < yi < 1. We define two embedded 
cylinders P, Q C as follows: 

P = {{x,y) : X e§>^,yo ^y i^yi} and Q = {(x, y) : xq < x ^ xi, y G S^}. 

We shall call P a 'horizontal' annulus and Q a 'vertical' annulus. We denote hy R = 
P U Q the manifold on which our linked-twist map will be defined and hy S = P D Q 
the 'intersection region'. See Figure [T7S1 

The set OPq = {{x,y) : x G S^,y = yo} denotes the 'lower' boundary of P, with the 
'upper' boundary dPi defined similarly. The 'left-hand' boundary of Q is denoted dQo 
and the 'right-hand' boundary denoted dQi. Again, these are defined similarly. Finally 
we denote dP = dPo U dPi and dQ = dQo U dQi. 



1 









Q 
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p 









xq Xi \ 

Figure 1.4: The manifold RcT^ (shaded). 

It is convenient to define the twist functions / and g from which our twist maps 
will be constructed on all of S^, as opposed to just on [yoiyi] and [xo,xi] respectively 
(as would most naturally fit in with the abstract definition above). Let / : S"*^ ^ be 
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given by 

{y- yo)/{yi - yo) if ?/ e bo, yi] 

fiy) = \ 

otherwise, 

and similarly g : S"*^ ^ by 



g{x) = < 



{x - xo)/{xi - xo) if X G [3;o,xi], 
otherwise. 



Both of these functions have the form (jl.2.ip illustrated in Figure ll.3f a) (recall that 
and 1 are identified in S^). They are differentiable for y G S^\{yoiyi} and x G 
S^\{a;o,xi} respectively. 

A horizontal twist map F : T'^ ^ is given by F{x,y) = {x + f{y),y) and it 
follows that F is continuous on and differentiable on T'^\dP. We remark that F is 
a homeomorphism of and is the identity map outside of P. We say that F is linear 
because of the piecewise linearity of /. 

Analogously we define a vertical twist map G : ^ by G{x,y) = {x,y + g{x)) 
and similar comments apply; in particular G = id outside of Q. 

A linear linked-twist map Hjj. : — > T'^ is given by the composition o F^ for 
positive integers j and k. We consider the restrictio n of Hj k to the invariant set R. 



Both twist maps pre serve the Lebesgue measure (see 



Katok and Hasselblatt (19951 ) or 



Sturman et al. (20061 )). so the composition Hj^k does also. We denote the Lebesgue 
measure on i? by ^. 

If we take xq = yo = and xi = yi = 1 and also j = k = 1 then Hj^k is precisely 
the cat map we have mentioned in Section [1.1. 1[ 



Finally, let us consider Hj j. as an abstract linked-twist map. ISturman et al. (20061 ) 
call the map co- twisting because jk is positive. We can take F = Ei o Ti o E^^ where 
Ti : Ci — > Ci is the linear twist map (defined by UfLTl} ) on Ci = x [yo,yi], and 
where 
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Ei{x,y) = {x,y) 



Similarly we have G = E2 o T2 o E2 where T2 : C2 ^ C2 is the linear twist map (again 
defined by (jl.2.ip ) on C2 = x [xq, xi] and where 

The Jacobians of Ei and E2 are given by 




respectively. The former has determinant 1 and the latter has determinant —1. The 
fact that the signs are opposite shows that Hj^k is co-twisting. 
Chapter [3] is devoted to proving the following: 

Theorem 1.3.1. If j and k are each at least 2, and one of them is at least 3, then 
the manifold C has an invariant Cantor set on which the linked-twist map Hj^k is 
topologically conjugate to a full shift on N = (j — 1)(A; — 1) symbols. 

1.3.2 Linked-twist maps in the plane 

Linked-twist maps in the plane have been studied by a number of authors. We provide 
the definition and state a theorem to be proven in Chapter SI In Section 12.3.21 we 
discuss the existing literature on planar linked-twist maps. In Sections 12.4.21 and 12.4.31 
we discuss some physical systems for which planar linked-twist maps provide a natural 
model. 

When dealing with linked-twist maps in the plane it will be convenient to denote 
§^ = [— vr, vr] where the opposite ends of the interval are identified. Let L be an annu- 
lus in the plane, centred at the origin and having inner and outer radii of vq and ri 
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respectively (where of course < ri), i.e. 

L = {{r, 9) : rQ ^ r ^ ri} 

where (r, 9) G M^j" x §^ are the usual polar coordinates. For convenience in what will 
follow, we assume that ri < vr. We observe that L is a cylinder as in our previous 
discussions. 

Define functions M± : M([ x ^ by 

M±{r, 9) = ±(r cos 6* — 1, r sin 9). 

The images M±{L) are annuli of the 'same size' in the plane, centred at (—1,0) and 
at (1,0) respectively. The annuli in the plane are expressed in Cartesian coordinates, 
which we denote by {u,v). Let A± = M±[L) denote these annuli. 

Under certain restrictions on ro,ri the annuli intersect in two distinct regions; this 
will be a necessary though not a sufficient condition for what follows and we will say 
more on the sizes of annuli later. We denote the intersection region in which the v 
coordinate is positive by S+ and the other by $]_. See Figure [L5l Let A = Aj^ U A_ 
and let S = U S_ . 



V 




\ , , , u 

-10 1 

Figure 1.5: The manifold A C (shaded). 
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Inverse functions ^ -.M? ^ M.^ x are given by 

A twist map A : L — > L is defined in polar coordinates: 

A(r, e) = {r,e + 27r(r - ro)/(ri - ro)). 

The twist function r 27r(r — rg) / (ri — rg) has derivative c = 27r/(ri — tq) and is affine; 
as before we abuse the notation shghtly and call it 'linear'. It has the form (11.2. ip 
illustrated in Figure [L3l fa). 

We define twist maps on A± as follows: let <I>,r : — > be given, respectively, 

by 

M+ o Ao M^^{u,v) ii{u,v)£A+ 
(u, v) otherwise. 



<I>(n, 



r{u,v) = < 



M_ o A-^ o MZ^{u, v) if {u, v) E A. 

(n, v) otherwise. 



A planar linked- twist map Q : A ^ Ais given by the composition G = Fo^. Figure [L6l 
illustrates its behaviour. 

We make some comments. First, we have defined as the composition of one twist 
map ^ and one twist map F, in contrast to our definition of a toral linked- twist map 
which was the composition of j 'horizontal' and k 'vertical' twists. We can of course 
define a more general planar linked-twist map F'^ o for j, k £ N. In fact all of the 
results we prove will go through with only trivial alterations; the cost however would 
be more cumbersome notation in a number of places. For this reason alone we take 
j = k = l. 

Second, the map G preserves the Lebesgue measure on A; see 



Woitkowski (19801 ). 
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Third, let us consider the planar linked-twist map as an abstract linked-twist map. 
The linked-twist map <I> restricted to is given by M-|_ o A o M^^ where M_|_ is the 
smooth embedding of cylinder L into the plane, and where A denotes the twist map. T 
restricted to A„ is given by o A^^ o AlZ^- 

It is convenient to express F in terms of A rather than A^"*^ (the latter not fitting our 
exacting definition of a twist map because the twist function has negative derivative). 
To this end we introduce a map B : L ^ L given simply by -B(r, 6) = (r, —6); it is easy 
to show that A^^ = B o Ko B^^. Our two embeddings are thus M+ and M_ o B and 
we will compare the signs of the determinants of their Jacobians in order to determine 
whether Q is co- or counter-twisting. We have 



DM±{r,0) =± 




and so determinants 

r (cos^ 9 + sin^ 6) , 

which clearly are both positive. It is easy to see that DB will have determinant — 1. 
Thus the embedding of L into preserves orientation whereas the embedding of L 
into A- reverses it; G is co-twisting. 

We illustrate the map's behaviour in Figure [TTUl 




Figure 1.6: One iteration of the planar linked-twist map. Part (a) shows some initial 
conditions in the form of a red horizontal line across the left-hand annulus A^. Part (b) 
shows the image of these points under the twist map <I> and part (c) shows the image 
under the linked-twist map = F o <I). 
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In Chapter m we will prove the following: 

Theorem 1.3.2. Let ro = 2 and ri = \/l ■ Then the planar linked-twist map @ : A —)■ A 
has the Bernoulli property, which is to say that it is isomorphic to a Bernoulli shift. 



Wojtkowski (19801 ) ■ in the particular case where the 



This verifies a conjecture of 
annuli are as stated. 

We comment that a weakness of our method is the need to be specific about the size 
of the annuli. We discuss this more in Chapter [6] where we are able to isolate which 
part of our proof would need to be improved upon to obtain a more general result and 
discuss our ideas for how this might be achieved. 

1.3.3 Linked-twist maps on the two-sphere 

In this section we introduce a linked-twist map on the sphere. One of the main features 
of interest is the construction of the pair of embeddings for which we make use of 
Jacobi's elliptic functions. We will state a theorem to be proven in Chapter [5j 

We review a small number of facts about Jacobi's ellip tic functions; we review severa l 
more in Section [53TJ For a comprehensi ve treatment see 

or alternatively see the excellent paper of lMever (200ll ). We plot the functions sn, cn, dn : 



Whittaker and Watson fl92m 



M — > [—1,-1] although we do not define them explicitly. Commonly each function de- 
pends upon a parameter k G (0, 1) also but we will always take k = \/2/2 so we omit 
this dependence. Let 

/'7r/2 

K{k)= / {1 - e sinH)'^^^ dt. 
Jo 



We comment that K = K{\/^/2) ~ 1.85. Functions sn and cn are periodic with period 
4:K whereas dn is periodic with period 2K; see Figure [TTTl 

It will be convenient to denote = [—2K, 2K] with the opposite ends identified. 
Let C = §^ X I where / = [—yo, yo] for some < yo < K. We define a second 'rotated' 
cylinder C = I' x where /' = [— xq, xq] for some < xq < i^. We can picture C, C 
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Figure 1.7: The Jacobi elliptic functions: sn(t) shown in red resembles a stretched sine 
function, whereas cn(t) in blue resembles a stretched cosine. The function dn(t) in green 
has no analogy amongst the standard trigonometric functions. 



as subsets of the two-torus = x as in Figure ll.8l fa) . 
Now let E : T"^ ^ S"^ C he given by! 

E{x,y) = (sn(j;)dn(?/), cn(x)cn(y), dn(x)sn(y)) . 

The restriction of E to either C or C' (though not their union) is a diffeomorphism 
between that cylinder and its image in §2 c (we prove this in Section 15.4. ip . Let 

= E{C), = E{C') and A = A+UA^. Our linked-twist map will be defined on A 
which is illustrated in Figure nTBT b) . (The fact that E is only injective when restricted 
to one of C or C' accounts for the fact that C Pi C" has one connected component but 

n Aj^ has two. We will say much more on this in Chapter [5j) 

As in our construction of an abstract linked-twist map we define a twist map on 
$ : Aj^ — > Aj^ to be the composition 

$ = EoToE~^, 
■^We are most grateful to Dr. Holger Waalkens for suggesting this map to us. 
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2K 



(a) 



-2K 





-2/0 



-2K 




Xq Xq 




C 



C 



2^ (b) 




Figure 1.8: It is simplest to define the manifold ^ C as the image of C U C" with 
respect to the function E defined above. Notice that our restrictions on the size of C, C 

lead to two distinct areas of intersection of the anmili Aj^ and A_ and four 'holes' which 
are not part of A. We have used dark shading for A and light shading for the 'reverse' 
side of A, as seen through the holes. 

where T : C — C is a linear twist map as defined previously. We extend $ to all of A 
by declaring it equal to the identity function on A- \Aj^ . 

One way to define a twist map T : A- A- would be to define a twist map on C', 
say T' : C' C , and define T = E oT' o E~^ . Instead we introduce a diffeomorphism 
N :C ^ C given by 



It is easy to check that N{C) = C. We define F : A- — A- as the composition 



and declare F equal to the identity function on A+\A^. The advantage of this definition 
is that it is trivial to show that the Jacobian matrix for N has negative determinant, 

and this leads immediately to the conclusion that if E is orientation-preserving then 
E o N must be orientation-reversing, and vice versa. Define the linked-twist map 




T = EoNoToN-^oE 



r-l 
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e = ro$, 



then Q is co-twisting. As with the hnked-twist map in the plane we could define a 
more general linked-twist map on the sphere as the composition F'^ o for j, k £ N. 
Our work would again require only trivial alterations but at a cost of more cumbersome 
notation. 

In Chapter Owe will prove the following. 

Theorem 1.3.3. The linked-twist map Q : A ^ A has the Bernoulli property, which is 
to say that it is isomorphic to a Bernoulli shift. 
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2 Literature Review 



The ergodic theory of hyperbolic systems is a significant branch of the dynamical sys- 
tems theory and the literature is appropriately rich and diverse. In this chapter we 
provide some key definitions and results on which our work builds, but in doing so we 
barely scratch the surface of all that is out there. 

The chapter is divided into four sections. In Section 12.11 we provide some basic 
definitions from ergodic theory, starting with the relatively weak property of ergodicity 
and building up to the strongest Bernoulli property. In Section 12.21 we review some 
key concepts and results from the hyperbolic theory. The systems we study in later 
chapters will all display non-uniform hyperbolicity and here we describe in detail what 
this means. 

In Section 12.31 we survey the literature pertaining to the linked-twist maps we de- 
scribed in the previous chapter. Finally in Section 12.41 we survey a few examples of 
applications for which linked-twist maps provide a natural model. The existence of such 
applications goes some way to explaining the recent resurgence in interest in linked-twist 
maps, as is perhaps best evidenced by the book of iSturman et al. (20061 ). 

2.1 Ergodic theory 

Ergodic theory is concerned with dynamical systems on measure spaces. It is typically 
highly non-trivial to prove that a given dynamical system has any of the ergodic prop- 
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erties we will present. However, the pay-off for doing so is a substantial amount of 
information about the behaviour of 'most' trajectories. 



2.1.1 Ergodicity and mixing 

All definitions and results in this section m ay be found in 



other standard reference for this material is 



Brin and Stuck C2OO2I I. An 



Katok and Hasselblatt (19951 ). 



Let (M, il, /, n) be a measure-preserving dynamical system. Here M is a set and will 
usually be furnished with some additional structure; typically we might require that M 
be a compact metric space, or a Riemannian manifold, il denotes a a-algebra of subsets 
of M, f a transformation of M into itself and fi a positive measure defined on il. 

Typically will be finite and so without loss of generality we may assume that it is 
a probability measure, i.e. /u(M) = 1. We will assume that / preserves /x, in the sense 
that for each set A E il we have fi{f~^{A)) = fi{A). 

Definition (Ergodicity). A dynamical system (M, il, /, /z) is said to be ergodic if when- 
ever ^4 G il has the property that f{A) = A, then either fi{A) = or /x(A) = 1. 

Ergodicity may be thought of as indecomposability, in the sense that two disjoint, 
non-trivial (i.e. positive measure) invariant sets are not possible. 
A stronger condition than ergodicity is the following: 

Definition (Strong mixing). A dynamical system (M, il, /, /x) is said to be strong mixing 
if for all sets A, B £ il one has 

hm /x(r"(A)ni?) =M^)Mi?). 

n— >oo 

The strong mixing (typically called just mixing) property implies ergodicity and 
can be thought of as points 'losing memory' of where they started. This is the kind of 
property we would like to prove for our dynamical systems. 
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In fact, it will be possible to prove a stronger property known as the Bernoulli 
property. The Bernoulli property is significantly more abstract than the other ergodic 
properties we have presented but the pay-off is substantial; Bernoulli systems behave, 
in a rigorous sense, as randomly as possible. 

2.1.2 The Bernoulli property 



Wiggins r2003l ). Let S = {1, 2, N} 



A good reference for the material in this section is 
be a collection of symbols, where N is an integer strictly greater than one. A (bi- 
infinite) symbol sequence has the form s = S-i, sq, si, ... where each Si G S. The 
space of all such symbol sequences is naturally thought of as the bi-infinite Cartesian 
product ■■■xSxSxSx---. We can define a metric on S^: if t = ...,t^i,tQ,ti, ... is 
another symbol sequence then let 



d{s,t)= ^ where 5i 

i=—oo 



if Si = ti 

1 otherwise. 



See 



Devanev (19861 ) for a proof that d is indeed a metric on S^. Intuitively points 
in T,^ are close if t heir sequences agree on a long central block. It is shown in 
Sturman et al. (20061 ) that the metric space ($]^,(i) is compact, totally disconnected 



and perfect (i.e. a Cantor set) and has the cardinality of the continuum. 

We now outline how one may define a measure on S^, following the approach of 



Arnold and Avez (19681 ). Let be the set of points in having j € 5" for the i*^ 
element in the symbol sequence. These sets generate a u-algebra of subsets of S^. We 
also define the cylinder sets 

h=l 

Define a normalised measure on S by insisting that for each j £ S we have ^ 



and J2j=i /^(j) = 1- The measure of a set is defined hy (Alj = and we extend 
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this measure to the cyhnder sets via the identity 



h=l 



It can be shown that ^ satisfies the axioms of a measure. 

The last part of our construction is a map a of into itself, known as a shift 
map. It is expressed most concisely by the relationship [o"(s)]i = Sj+i, although perhaps 
intuitively it is preferable to insert a period at some point in the symbol sequence 
(written without commas), and look at where that period occurs in the symbol sequence 
of the image: 

S = • • • S-2S-I.S0S1S2 • • • , Cr(s) = • • • S-2S^lSo.SlS2 • • • . 



If the domain is all o f then a is often called a full shift on N symbols. In this case 
Wiggins (2OO3I )) that a is a homeomorphism of S^, that it has 



it can be shown (see 



a countable infinity of periodic orbits including orbits of all p eriods, an uncountab le 



infinity of non-periodic orbits, a dense orbit and moreover (see lSturman et al. (20061 )) 



that a preserves the measure fi constructed previously, with respect to which it is mixing. 

Definition (Bernoulli property). A dynamical system (M, it, /, z^) is said have the 
Bernoulli property if it is (metrically) isomorphic to a Bernoulli shift. More formally, 
we require that the following diagram commutes 

a 



M 



4> 



M 



where cp : {T,^ , fi) ^ (M, is an isomorphism. 



A map having the Bernoulli property automatically has all of the properties of 
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the full shift on symbols, given above. An example of a dynamical system hav- 
ing the Bernoulli pr operty is the well-known baker's map of the unit square. See 



Sturman et al. (20061 ) for further details of the map and a proof of this result. In this 



example the isomorphism can be constructed explicitly; this is aided by discontinuities 
in the map, in contrast to the maps we will consider. 



2.2 Hyperbolicity 

Hyperbolicity is an important part of the dynamical systems theory and the focus of a 
great deal of active research. Knowing that a certain dynamical system has hyperbolic 
structure gives us access to a number of results and techniques for demonstrating ergodic 
properties. All of the dynamical systems we consider in this thesis display hyperbolicity; 
here we outline some of the key definitions and results we will repeatedly rely upon. 

2.2.1 Uniform hyperbolicity 

Let (M,it, /, /i) be a dynamical system. Throughout this section we will assume that 
M has some structure beyond being just a set, likewise / some smoothness properties. 
In particular we assume that M is a compact, smooth (C°°) n-dimensional Riemannian 
manifold and / a smooth (C°°) diffeomorphism of M (i.e. a differentiable map with 
differentiable inverse; we have defined these terms in the case n = 2 in Section [1.2. ip . 
In section [2.2.3l we will discuss what happens when we relax these conditions somewhat, 
but for now let us keep the exposition as clean as possible. Denote by Df^ the Jacobian 
matrix of / evaluated at x £ M. 

Definition (Hyperbolic fixed point). Let x £ M be a fixed point of f , i.e. f{x) = x. 
Then x is said to be hyperbolic if none of the eigenvalues of Dfx have magnitude one. 

In some neighbourhood U of a hyperbolic fixed point x £ M, the dynamics of / 
will closely resemble the behaviour of the linearised system. To be precise, there is a 
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corresponding neighbourhood V of the origin and a homeomorphism h : V ^ U such 
that f(h(i/)) = hiPfxU) for ah y £ V. This is the well-known Hartman-Grobman 



theorem; see 



Robinson (Wm . 



More generally we will define hyperbolicity on a set, rather than at a fixed point. 



The simple st case is where the 
taken from lBrin and Stuck (20021 ). 



lyperbolicity is uniform. The following definition is 



Definition (Uniformly hyperbolic set, Anosov diffeomorphism) . Let {M,ii,f,^) be a 
measure-preserving dynamical system. Let U <Z M he a non-empty open subset and 
f : U ^ f{U) C M a smooth diffeomorphism. A compact, f -invariant set A C U is 
said to be uniformly hyperbolic if there exist constants c > and < A < 1, and if there 
is a continuous splitting of the tangent space T^M = E^(x) (B E^(x) at each x G A such 
that 

Df,E%x) = E^fix)) and Df,E^{x) = E^{f{x)), (2.2.1) 



||L»/>,|| < cA"||i;,|| for £ E'ix) (2.2.2) 

\\Df-^vJ ^ cX^VuW for vu G E^{x). (2.2.3) 

If A = M then f is called an Anosov diffeomorphism. 

Condition (j2.2.ip says that stable and unstable directions should be invariant un- 
der the differential map Dfx, whereas conditions (j2.2.2|) and (j2.2.3|) give estimates on 
the contraction of stable subspaces under forward iteration, and of unstable subspaces 
under backward iteration respectively. We call this hyperbolicity 'uniform' because the 
constants c and A are independent of the point x G M. Notice that this is precisely the 
situation we found when analysing the cat map in Section II. 1.11 
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It turns out that few dynamical systems are uniformly hyperbolic. In the next 
section we discuss a generalisation. 



2.2.2 Pesin theory 

We begin by defining what we mean by non-unifo rmly hyperbolic, then go on to de- 



scribe some celebrated results due to 



Pesin (19771 ) which have hugely influenced the 



stu dy of such systems oy er the past three decades or so. 



Barreira and Pesin (2002 ) 



and ISturman et al. (20061 ) both give good accounts of the material in this section; our 
definitions are taken from these sources. 

We will not use the results from this section as they do not apply to the linked-twist 
maps we study, which are not diffeomorphism s. However they pr ovide a natural bridge 



Katok et al. (19861 ) that we will describe 



from uniform hyperbolicity to the theorem of 
next and use extensively thereafter. 

Definition (Non-uniformly hyperbolic). The measure-preserving dynamical system (M, il, 
is said to be non-uniformly (completely) hyperbolic if there exist measurable functions 
< A_(x) < 1 < A-i-(x) and e{x) such that e(x) = e(/(x)) (i.e. e is invariant along 
trajectories), and if there is a splitting of the tangent space TxM = E^[x) © i?"(x) for 
each X, and finally a function c{x) so that, for each k £ Z and n > we have 



Df^E'{x) = E'ifix)) and Df^E'^ix) = ^"(/'^(x)) 



7^U( fk/ 



(2.2.4) 



\\Df>s\\ ^ c {^f\x)) \l[x)\\v,\\ for V, € E'{x) (2.2.5) 
Df^VuW ^ (^f\x)) Xl{x)\\vu\\ for G E^{x), (2.2.6) 



29 



Z{E'{x),E^{x))^c"\x), (2.2.7) 

^ c(x)e"(^)l^l. (2.2.8) 

Conditions (|2.2.4p . (j2.2.5p and (j2.2.6p are analogous to the conditions we impose on 
uniformly hyperbolic systems, though the replacement of independent constants with 
functions of x means that they are less restrictive. Condition (j2.2.7p says that the 
stable and unstable directions are transversal. The final condition (j2.2.8p is perhaps 
a little more subtle and deals with the rate at which our contraction or expansion 
estimates in conditions (j2.2.5p and ()2.2.6p deteriorate along a trajectory. It says that this 
deterioration is sub-exponential and is thus dominated by the exponential contraction 
or expansion. 

The most important tool for analysing non-uniformly hyperbolic systems is the 
Lyapunov exponent. 

Definition (Lyapunov exponent). For a dynamical system (M, il, /, /x), the Lyapunov 
exponent x^(^) ^) point x £ M and in the direction v S TxM is given by 



X^{x,v)= lirn -log||D/>||, 

n— >±oo n 



whenever this limit exists, where \\ ■ \\ denotes the standard Euclidean norm in M"'. 

The importance of Lyapunov exponents is illustrated by the fact that non-uniformly 
hyperbolic systems are commonly known as systems with non-zero Lyapunov exponents. 
The result to which this epithet alludes is the following. 



Theorem 2.2.1 (jPesin (19771 )). A dynamical system (M, it, /, /i) is non-uniformly (com- 
pletely) hyperbolic if for almost every x €z M the Lyapunov exponent xi^i v) is non-zero 
for every non-zero v £ TxM . 
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Pesin derived two further results which lay the foundations for our results on linked- 
twist maps. The first is the famous stable man ifold theorem. We need the following 



definition, taken from 



Barreira and Pesin (2002 ). 



Definition (Local invariant manifolds). Let i?*(0,e) be the open e -neighbourhood of the 
origin in E^{x), and similarly B^{0,e). A local stable manifold of x has the form 



7^(x) = exp,.{(x,V'^(:E)) : x G B-'{0,e)} 



(2.2.9) 



for some £ > 0, where : B^{0,£) — > E'^{x) is a smooth map satisfying ip'^{^) = and 
Dtl'^{0) = 0. The trajectories of x and y G 'J^{x) approach each other at exponential 
rate as n ^ +oo. Transposing u and s in the above and considering n ^ —oo yields 
the description of local unstable manifolds; we omit further details. 



Theorem 2.2.2 (|Pesin (19771 )). If f : M ^ M is of class C"^"'"" and is non-uniformly 
hyperbolic, then for almost every x £ AI there exists a local stable manifold 'y^{x) with 
the properties that x £ 7*(x), Tx^^{x) = E^{x) and if y £ "y^ix) and n ^ then 



p(r(a;),r(y))^T(x)A"e-p(x,y), 

where p is the distance in M induced by the Riemannian metric and T : M ^ (0, oo) 
is a Borel function satisfying 

r(/'"(x)) ^ r(a;)e^°"H^ m G Z. 



The other result that will be crucial to our work shows that the manifold M has an 
ergodic partition, the definition of which is given in the following statement. 



Theorem 2.2.3 (jPesin (19771 1). If f : M ^ M is of class C^"*"" and is non-uniformly 



hyperbolic then M is either a finite or countably infinite union of disjoint measurable 
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sets AIq, Ml, ... such that /i(Mo) = and fi{Mi) > for all other subsets, each Mj is 
f -invariant (i.e. f{Mi) = Mi) and the restriction of f to any Mi is ergodic. 

This completes our very brief exposition of Pesin theory. As we have mentioned, 
in order to use Pes in-type results as above we need to appeal to more general work of 



Katok et al. (19861 ) which we review in the following section. 



2.2.3 Smooth maps with singularities 



This section surveys the resul ts of 
that found in the appendix of 



Katok et al. f 19861 ). our exposition following closely 
Przvtvcki (198j ). In stating the results it is necessary to 



introduce several nested full-measure sets. We tabulate these (Table 2.1) in the hope 
that it helps the reader through the construction. 

Let X be a complete metric space with a metric p. Let M C X be an open 
subset which is also a Riemannian manifold, with the Riemannian metric inducing p\m- 
Assume there is some r > such that for each x G M the exponential map exp^,, 
restricted to the ball 



B{x) = B{x,min{r,disip{x,X\M))), 



is injective. 

Let (M, il, /, /i) be a measure-preserving dynamical system as before, but with the 
difference that we define / only on an open set N C M into M. The measure p, is an 
/-invariant probability measure on M and we require that, on A^, the function / is 
and injective. Finally we denote sing(/) = AI\N. 

We call / a smooth map with singularities. This completes the definition of the map 
itself. We now describe two conditions of a technical nature that place restrictions on 
the nature of the points at which Df is undefined, and on the growth of D^f near to 
these points. 
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Set 



X 

M 
N 
J 
E 



Description 



Complete metric space 

n-dimensional Riemannian manifold 

Open set on which / is defined 

Intersection of all images and pre-images of N 

Set on which Lyapunov exponents exist 



Table 2.1: Full measm'e subsets of the complete metric space X. As we have listed the 
sets, each contains those below it. 



In keeping with the notation of 



Przvtvcki Mm and 



Sturman et al. (20061 ) we say 



that / satisfies the condition (KSl) if and only if there are positive constants a and ci 
so that for every e > we have 



/i(i?(sing(/),e)) ^Cle^ 



(2.2.10) 



where i?(sing(/), e) means the open e-neighbourhood, with respect to p, of the set 
sing(/). 

We say that / satisfies the condition (KS2) if and only if there are positive constants 
h and C2 so that for every x £ N we have 



\D^f{x)\\^C2p{x,smg{f)) 



-b 



(2.2.11) 



where ||D /(x)|| denotes the supremum of \\D {exp'^ of o expy)||, taken over those x 
for which x £ B{y) and /(x) G B{z). 

Informally, the two conditions state that 'most' points have neighbourhoods free 
from singularities and that the second derivative does not get large 'too quickly' (e.g. 
exponentially) as we approach these singularities. 

If / satisfies condition (KSl) then /i(sing(/)) = and so p{N) = 1. Let J = 
n^_oo/"(-^)- '^^^ /-invariance of p implies that p{J) = 1 also ( we prove this in 



Chapter [5]). The Multiplicative Ergodic Theorem (originally proven by 



Oseledec (1968 ) 
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but see also IChernov and Markarian (20031 ) and iBarreira and Pesin (20021 )) holds for 
smooth maps with singularities satisfying the above conditions. We have 

Theorem 2.2.4 (jChernov and Markarian (20031 )). Suppose that 



[ log+ p/IUpd^ < oo and [ log+\\Df-^\\apd^i<oo, (2.2.12) 
Jx Jx 

where fog^(-) = max{fog(-), 0} and || • ||op denotes the operator norm induced by p. Then 
there is an f -invariant set E C J, /u(-E') = 1, so that for every x G E and every non-zero 
V G T^X , the Lyapunov exponents x^{x,v) exist. 

The following theorem provides the fr amework for our wor k in Chapters HI and El 
in addition to the main reference, see also 



r 

Przvtvcki fl983 ) and .Sturman et al. (20061 ). 



The theorem contains the definition of an ergodic partition to which we will refer several 
times, and includes a more complete description than was given in Theorem 12.2.31 



Theorem 2.2.5 ( Katok et al. (19861 )). Let f be a smooth (C^ at least) map with sin- 



gularities, defined a.e. on smooth manifold M as above. 

(a) Suppose f satisfies the conditions (KSl) and (KS2) and the hypothesis of The- 
orem \2^2.4\ above. Then for a.e. x S M and for all non-zero tangent vectors 
V G TxM , the Lyapunov exponents x^i^^v) exist. Corresponding to any posi- 
tive (respectively, negative) Lyapunov exponents, there exist local unstable (stable) 
manifolds 7"^*) (x) of the form we have described. 

(b) // additionally for a.e. x £ M and for all non-zero v £ TxM we have x{x,v) 7^ 
then M decomposes into an (at most) countable family of positive measure, 
f -invariant, pairwise- disjoint sets Mi on which the restriction of f is ergodic. 
Furthermore each set Mi has the form Mi = |J"f!| M^ where, for each j, 

is Bernoulli. Such a system will be said to have an ergodic partition. 

(c) // additionally for a.e. x,y £ Al there exist integers m,n such that 
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r(7"(x))nr"(7^(2/))/0 



(2.2.13) 



(we say that f satisfies the manifold intersection property), then in the decompo- 
sition of M there is just one positive-measure set, i.e. Mi = M for all i, and so f 
is ergodic. 



(d) // additionally, for a.e. x,y £ M, the condition ^2.2.13\) is satisfied for each pair 
of sufficiently large integers m,n (we say that f satisfies the repeated manifold 
intersection property), then f has the Bernoulli property. 



2.2.4 The Sinai-Liverani-Wojtkowski approach to proving ergodicity 



We discuss some work of 



Liverani and Woitkowski (19951 ) based on work of lSinai (19701 ) 



le class of sy s 



NicoUl996a|)) 



in which the authors estabhsh criteria for certain maps to be ergodic. Tl 
terns to which their results apply is broad and includes some (according to[ 
of the linked-twist maps we will study, if not all of them. There are a large number 
of technical hypotheses in the statement of their theorem and we do not intend to use 
their results, so we limit our exposition to an informal discussion. 

In part we do not appeal to their result because the aforementioned technical consid- 
erations mean that this would be far from trivial (although this alone is not justification, 
as the same can be leveled at our approach!). However, when we come to make our clos- 
ing remarks and discuss our ideas for future works we will have good cause to conjecture 
that our approach offers benefits t hat theirs do es not. 



The significance of the works of 



Pesin (19771 ) and later of 



Katok et al. (1989 ) is that 



they extended the class of sys tems for wh ich an erg odic partition (a local property ) 



can be established. Conversely 



Sinai (1970 ) and later 



Liverani and Woitkowski (19951 ) 



extend the class of systems for which ergodicity (a global property) can be established. 

Consider the condition (j2.2.13p given in Theorem 12.2.51 which is the bridge be- 
tween local and global properties in that theorem. The condition is sufficient because 
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it allows one to conclude that any integrable observable on M (that is, any function 
g £ (M, M)) that is /-invariant, is constant //-almost eve rywhere. It can be show n that 



this condition is equivalent to ergodicity; see for example 



idea dates back to 



Brin and Stuck (2002 ). This 



Hopf (19391 ) and such a construction is s ometimes called a Hopf chain; 



Liverani and Woitkowski (19951 ). 



for more details we recommend the introductory sections of|] 

Proving that (j2.2.13p is satisfied, in general, will require some specific knowledge 
of the nature of local invariant manifolds. If we are able to conclude that the length 
of f^{l^{x)) diverges with m and that the length of / "(7'^(y)) diverges with n, and 
moreover if we can give a 'useful' characterisation of their respective orientations, then 
we might hope to satisfy the condition. Either, or both, of these may be far from trivial 
within a generic non-uniformly hyperbolic system. 

As we have seen, such systems do not have uniform growth rates for local invariant 
manifolds, nor indeed uniform lower bounds on the original sizes of those manifolds. 
Thus in general one needs a more sophisticated approach, and this is precisely where 
the Sinai-Liverani-Wojtkowski approach comes in. 

The Sinai-Liverani-Wojtkowski method (as we call it) works by constructing a con- 
nected 'chain' of local invariant manifolds with one end at x and the other at y. So far 
this is just Hopf 's method, but rather than relying on growth and orientation to deduce 
this connection, their method is to very carefully partition M using small overlapping 
'squares' whose sides are parallel to the stable and unstable directions respectively. 

Their arguments relate the width of a chosen partition to the conditional measure 
of those points within a given square whose local invariant manifolds completely cross 
that square. In this manner they are able to conclude that if a point x G M satisfies 
a certain local condition, then there is an open set containing x that is itself contained 
within a single ergodic component. It is an easy corrolary that if /x-a.e. x G M satisfies 
this condition then the map is ergodic. 

As we have briefiy argued, the real achievement of these methods is to deduce 
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ergodicity in systems where either the growth or orientation (or both) of local invariant 
manifolds are not well-behaved, in some sense. The linked-twist maps for which we 
prove strong ergodic properties do not fall into this category; we think in particular 
of the linked-twist map in the plane: that local invariant manifolds grow arbitrarily 
long is known and mentioned elsewhere; that the orientations of these manifolds can be 
characterised in a useful manner is the cornerstone of our proof. 



2.3 Linked-twist maps 

The linked-twist map literature spans almost three decades and includes results de- 
scribing certain ergodic properties of the toral and the planar linked-twist maps we 
have mentioned. In Section 12.3.11 we review the relevant results for toral maps and in 
Section 12.3.21 we do the same for t he planar maps. For a comprehensive overview of this 



literature the reader is directed to 



Sturman et al. (20061 ). In Section [2.3.31 we list some 



other explicitly defined maps for which strong ergodic properties have been established, 
in the hope that this will help to put the results for linked-twist maps into context. 

2.3.1 Linked-twist maps on the two-torus 

Let H = o : R ^ R, with C T^, be a toral hnked-twist map as defined in 
Section ll.3.1[ Recall that the product jk is positive if and only if a toral linked-twist 
map Hj k is co-twisting. The following theorems describe the ergodic properties of these 
maps. 



Theorem 2.3.1 (jBurton and Easton (198d )). If jk > and H is composed of smooth 



twists, then H has an ergodic partition. 

The smooth twists are as depicted in Figure [L3l fbl. Furthermore the authors sketch 
a geometrical argument with which this may be extended to the Bernoulli property. 
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Theorem 2.3.2 (jPevanev (19801 )) ■ If jk > then periodic and homoclinic points of H 
are dense, and H is topologically mixing. 

Devaney's result is more topological in nature; in fact this theorem was motivated by 
the similarities between toral linked- twist maps and the cat map. We observe also that 
Devaney does not need to put restrictions on the nature of the twist functions (aside, 
of course, from those conditions mentioned in Section 11.2.21 which all twist functions 
must satisfy). 



Theorem 2.3.3 (jWoitkowski (1980 )). If jk > and H is composed of linear twist 
maps, then H is Bernoulli. Alternatively, if jk < —4 then H has an ergodic partition. 

Wojtkowski's result on the torus is the main source of inspiration for our proof of 
the Bernoulli property in the pla ne. In fact he proved only the i^T-property; it follows 



from lChernov and Haskell (19961 ) that the system is necessarily Bernoulli. Wojtkowski 
also considers planar linked-twist maps in the same paper; we mention this in the next 
section. 

In order to state the next result we briefly describe what is meant by the strength 
of a twist. Consider a 'horizontal' twist map : P ^ P as defined in Section [1.3.11 
i.e. F^{x,y) = {x + jf{y),y) where / : / ^ S-*^ is a twist function. The strength of the 
twist map is defined to be 



SF = sgn(j) inf \jf'{y)\. 



Denote / = [«o,ii]- In the common case where / is affine then sp = j/{ii — io)- The 
strength of a 'vertical' twist map is defined analogously. 



Theorem 2.3.4 (|Przvtvcki (;i983l X^. If spsa < -Co ~ -17.24445 and \j\,\k\ ^ 2 then 
H is Bernoulli. 

Przytycki constructs an intricate argument allowing him to prove this result for 
counter-twisting maps. We have mentioned before that these are more difficult to study 
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than their co-twisting counterparts; a comparison of the criteria in this theorem and 
the previous one exemphfies the situation. 

Two other results for toral hnked- twist maps are known to us, though they take us 
a httle furt her afie ld so we do not state them precisely. Both are due to Nicol. In the 



first paper (|l996al ) he constructs a linked-twist map which has the Bernoulli property, 



despite having local invari ant ma nifolds and positive Lyapunov exponents only on a null 



set. In the second paper (|l996bl ) he considers a Bernoulli linked-twist map of infinite 
entropy having smooth local invariant manifolds and positive Lyapunov exponents a.e. 
with some discontinuities. He shows that the map is stochastically stable. 



2.3.2 Linked-twist maps in the plane 



Let e 



: A ^ A, with A C 



be a planar linked-twist map as defined 



in Secti on 11.3.21 The study of these maps was motivated by a number of authors. 
Thurston (198S ) encou ntered linked -twist maps such as these in his study of diffeomor- 



phisms of surfaces and 



Braun (198ll ) showed that simi lar maps arise as an approximate 



model of the global flow for the Stormer problem. iBowen (19781 ) showed that cer- 



tain linked-twist maps on such a manifold have positive topological entropy, and asked 
whether they possessed any ergodic properties. The following results describe what is 
known. As before jk > is the co-twisting case. 



Theorem 2.3.5 (jPevanev (19781 )). For non-zero j, k there is an invariant zero-measure 



Cantor set on which the map @ is topologically conjugate to a subshift of finite type. 

Devaney's paper motivates our construction of a similar invariant set for a toral 
linked-twist map. 



Theorem 2.3.6 (jWojtkowski (198d )). If jk > and the twists are sufficiently strong, 
then has an ergodic partition. Alternatively if jk < and a different (stronger) twist 
condition is satisfied, then also Q has an ergodic partition. 
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We give further details in Section [4.11 As we have mentioned, Wojtkowski's work is 
th e main source of inspiration for our proof. We mention also some unpublished notes 



of 
In 



Przvtvcki fl98ll ): the planar linked-twist maps are amongst those that he discusses. 



Przvtvcki (19861 ) the author considers again this large class of maps and shows that 



under certain conditions periodic saddles and homoclinic points are dense, and that the 
maps are topologically transitive. 



2.3.3 Other maps with strong ergodic properties 

Finally we describe some other systems for which strong ergodic properties have been 
established. The list, although not exhaustive, contains the majority of examples of 
which we are aware. As such it demonstrates the relative scarcity of such results and 
serves to underline the significance of our constructions. 

The main classes of examples and specific ex a.mples of which we are aware are: 



geodesic flow s on manifolds with negativ e curvature (lAnosov and Sinai (19671 ) , 



Burns (19831 )): gases of hard spheres 



pseudo-Anosov diffeomorphisms (|Anosov and Sinai (19671 ) 



Ballmann and Brin (19821 ). 



Kramli et al. f 19891)'): svmplecti c Anosov and 



Gerber a985l'l. 



geodesic flow s on surfaces with special metri cs and potentials (jPonnav (19881 ). 



Knauf ggSTl'l): svstems like Wojtkowski's (ll99C 



MacKav (2006l)): 



Burns and Gerber (19891 ). 



Barnes and Koss (20001 )): and the Belykh map 



; certain rati onal maps of the sphere 



Sataevjl999|)) 



Katok (19791 ) who showed that Bernoulli 



Moreover we mention the important work of|] 
diffeomorphisms may occur on any surface. Beginning with a hyperbolic toral automor- 
phism similar to that which we have encountered (and which is uniformly hyperbolic), 
Katok 'slows down' trajectories in a neighbourhood of the origin. The manner in which 
this slowing down is accomplished is a little sophisticated a.nd we do not intend to pr o- 



Barreira and Pesin (20021 ). 



vide the details here; an excellent account can be found in 

One of the consequences is that a Bernoulli map, derived from an Anosov diffeomor- 
phism, can be shown to exist on . This is particularly interesting given the well-known 
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results of lHirsch (197ll ) and lShiraiwa (19731 ) that cannot support an Anosov diffeo- 
morphism. 



2.4 Applications of the linked- twist map theory 



In recent years the study of hnked-twist maps has taken on a new significance owing to 
developments in our understanding of the mechanisms underlying good mixing of fluids. 
Ottino (19891 ) has shown that the single most important feature to incorporate in the 
design of any fluid mixing device is the 'crossing of streamlines', by which we mean that 



flow occurs periodically in two transversal directions. That 
suitable paradigm for this design process was hi ghlighted in 



and has been discussed at much greater length in 



inked-twist maps provide a 



Ottino and Wiggins ('20041') 



Wiggins and Ottino (20041 ) and 



Sturman et al. f2006l ^ 



2.4.1 DNA microarrays 

One important example of physical systems that may be analysed within the linked-twist 
map framework are cer tain models of DNA m icroarrays. In this section we summarise 



some ideas contained in 



Hertzsch et al. (20071 ): see that paper and the references therein 



for further details. DNA microarrays have been used widely in biochemical analysis for 
a number of years. Amongst their uses are gene discovery and mapping, gene regulation 
studies, disease diagnosis and drug discovery and toxicology. 

A DNA microarray consists of DNA strands ('probes') fixed to a surface such as 
glass or silicon. This array is placed in a hybridization chamber containing a solution of 
DNA or mRNA (the 'target'). Hybridization occurs when the target strand combines 
with a complementary probe strand, as governed by base-pairing rules. 

Hybridization is most efficient when each target strand can move throughout the 
solution and encounter every probe. Two processes lead to this: diffusion and advection. 
The former cannot be relied upon to produce the desired result in reasonable time 
because the typical situation involves low Reynolds number and thus no turbulence. 
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Advection in such devices has consequently been the focus of a great deal of research 
into how one might induce good mixing. 

It is now well established that chaotic advection provides a source of efficient mix- 
ing in many fluid problems, in particular those on the 'microfluidic' scale of the DNA 
microarrays. Two designs for such devices, both relying up on cyclic removal and rein- 



jection of fluid into the mixing chamber, are detailed in 



McQuain et al. (20041 ) and 



Raynal et al. (20041 ). Typically two different source-sink pairs are used. 

Two factors which have a great impact on the efficiency of such mixers are the 
locations at which fluid should be removed and reinjected, and the time for which such 
a source-sink pair should be active. Here linked-twist map theory can help to inform 
the design. The motion of fluid in such a device can bear striking resemblances to the 
motion of a point in the domain of a toral linked-twist map. Analysing these mixing 
devices in this manner has lead to the proposal of new mixing protocols. 

2.4.2 Channel-type micromixers 

There are many areas of applications where the linked-twist maps which most naturally 
act as models are defined on surfaces other than the two-torus. One such example are 
channel- type microm ixers. In this section we briefly summarise some ideas presented in 



Stroock et al. (20021 ): see that paper and the references therein for further details. 

Mixing of the fluid flowing through a microchannel is highly desirable in a number 
of situations, including the homogenisation of solutions of reagents in chemical reac- 
tions and in the con trol of dispersion of material along the direction of Poiseuille flows. 



Stroock et al. (20021 ) argue that, at low Reynolds number and using a 'simple' chan- 
nel (i.e. one that is straight and has smooth walls) any mixing is a consequence of 
diffusion only. Moreover they conclude that the rate at which this happens, even in a 
microchannel, is slow compared with convection along the channel. To reduce the length 
of channel required for mixing to occur one needs to introduce transversal components 
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of the flow which stretch and fold volumes of fluid over the cross section of the channel, 
thus reducing the distance over which diffusion must act. 

Such transversal flows may be generated by placing ridges on the floor of the channel, 
at an oblique angle to the flow. The ridges present anisotropic resistance to viscous flows 
resulting in transversal flow, which then circulates back across the top of the channel. 
Consequently flow along the channel becomes helical. The helical motion of the flow is 
a motion that can be approximated by certain planar linked-twist maps. 



2.4.3 Other examples 

The existence of the above examples illustrates that the linked-twist map framework can 
be a useful tool in the development of models for certain mixing devices. In the case of 

lere are nu merous other examples we could have mentioned. 



planar linked-twist maps t 
the blinking vortex flow of lAref (19841 ) being a prime case in point. Here a pair of point 
vortices in an unbounded inv iscid fluid are alterna tely switched on. Further work on 
this system was conducted bv lKhakhar et al. (1980 ). 



The work has applications t o the study of tidal flow close to a headland .jutting 



Signeh and Gever ri99lh . 



Signell and Butman (1992 ) 



out into the sea. Fo r details see 
and ISamelson (19941 ). In the latter reference a kinematic model to study mixing and 



transport by eddies is developed. For details on how linked-t wist maps may b e used as 



a paradigm for studying such systems we direct the reader to 



Wiggins (19991). 



Yet more examples are gi ven by the electroos motic stirrers of 



the cavity flows introd uced by 



and 



Qian and Bau f2002l'l: 



Chien et al. (19861 ) a nd developed further by 



Jana et al. (19941 ): and the egg-beater flows of 



Leong and Ottino (198% 



Franiione and Ottino (1992 ) 



Linked-twist maps on the two-sphere do not lend themselves to applications as read- 
ily as their counterparts in the plane, but we are able t o extract an example from the 



field of quantum ergodicity. 



Marklof and O'Keefe (20051 ) have shown that the quantum 



eigenstates of linked-twist maps defined on the two-torus are equidistributed. This 
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result uses th e fact that the corresponding classical linked-twist maps are ergodic. 



O'Keefe (20051 ) demonstrates that an analogous result can be shown to hold on the 
two-sphere, if one is able to show that the corresponding classical maps are ergodic. We 
mention another possible application for which linked-twist maps on the sphere might 
be a useful analytical tool in Chapter [6l 

Finally we mention an application to granular mixing. There are numerous situ- 
ations in the pharmaceutical, food, chemical, ceramic, metallurgical and construction 
industries where an understanding of the behaviour of granular media is crucial. How- 
ever the literature dedicated to the mixing of these materials is not nearly as developed 
as its counterpart for fluids. Recent work has shown that this is yet another example of 
a p hysical process which may be modelled using linked-twist maps. For further details 



see 



Sturman et al. (2008 ) and the references therein. 
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3 A horseshoe in a toral 
hnked-twist map 



This chapter is motivated by work of 



Devanev (19781 ) in which the author estabhshes 



the existence of a topological horseshoe in a planar linked-twist map. We construct, 
using similar ideas, the counterpart for a toral linked-twist map. To our knowledge this 
is not to be found in the literature. 

Let J, k be positive integers such that N = (j — l){k — 1) ^2 (i.e. each of j, k is at 
least 2 and one is at least 3), and let oF^ : i? — > i? be a toral linked-twist map, as in 
Section [1.3.11 We show that there exists a zero-measure, compact, invariant set within 
i? C (a 'horseshoe') on which the dynamics are topologically conjugate to a full shift 
on symbols. (By contrast, Devaney's construction yields a conjugacy with a suhshift 
instead; we discuss this further in Chapter [6l) 

In Sectio n 13.11 we develop some notation with which to state a theorem due to 



Moser (19731 ). which provides sufficient critera for the existence of the horseshoe. In 
Section 13.21 we show that toral linked-twist maps as above satisfy these criteria; this 
entails a detailed geometrical construction. Finally in Section 13.31 we provide some 
extra analysis to show that the map restricted to the horseshoe is uniformly hyperbolic. 

Before we begin we remark on our notation. It will be natural in this chapter to 
reserve the letter H for horizontal strips, which we define below. To avoid confusion, 
throughout this chapter, the linked-twist map will always be denoted by o (and 
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not by Hj is), whereas H,Hi etc. denote horizontal strips. 

3.1 The Conley-Moser conditions 

Here we describe sufficient criteria for a two-dimensional invertible map to possess an 
invariant Cantor set on which the aforementioned conjugacy exists. These are commonl y 



known as the Conley-Moser conditions, having first been introduce d iniMoser (lOTa ). 
A lucid and comprehensive account of this material may be found in lWiggins (20031 ). 

3.1.1 Horizontal and vertical curves and strips 

We begin with some definitions. Recall that a real-valued function / defined on con- 
nected domain C M is Lipschitz continuous if and only if there exists a constant c > 
and for every pair a,b £ D we have 

|/(a)-/(6)Kc|a-6|. 

We will say that such an / is c-Lipschitz. We use this notation to define curves in 
S = P n Q C R. Recall that S = [xo,xi] x [uq, yi]. 

Definition (m/j-horizontal and m^,-vertical curves). An rrih-horizontal curve is the 
graph of an rrih-Lipschitz function h : [xq,xi] [yo,yi]. An rriy-vertical curve is the 
graph of an m^-Lipschitz function v : [i/q, yi] [xq,xi]. 

We use such curves to form the boundaries of strips as follows. 

Definition (m/j-horizontal and mt,-vertical strips). Given two non-intersecting ruh- 
horizontal curves of functions h\ and h2, with hi(x) < h2{x) for each x G [xo,xi], an 
mfi-horizontal strip is the set 

H = {{x,y) £ S : X £ [xo,xi],y £ [/ii(x), /i2(x)]} . 
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The mfi-horizontal curves are then referred to as the horizontal boundaries of H . The 
vertical boundaries of H are contained within the lines x = xq and x = xi. 

Similarly given two non-intersecting my -vertical curves of functions vi and V2, with 
vi{y) < V2{y) for each y £ [yo,yi], an ruy-vertical strip is the set 



The my-vertical curves are then referred to as the vertical boundaries ofV. The hori- 
zontal boundaries of V are contained within the lines y = yo o,nd y = yi. 



Figure 3.1: The region S C R showing horizontal curves of y = hi{x) and y = h2{x), 
vertical curves of x = vi{y) and x = V2{y), a horizontal strip H and a vertical strip V. 

We illustrate some curves and strips in Figure 13.11 Define the width of a strip as 
follows: 

Definition (Width of m/i-horizontal and m„-vertical strips). Let H be an nih-horizontal 
strip as above. Its width is given by 



V = {{x,y) e S :y e [yo,yi],x g [vi{y),v2{y)]} . 




x = vi{y) 



x = v2{y) 



d{H) 



max \h2{x) — hi{x)\. 
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Similarly let V be an m^-vertical strip as above. Its width is given by 

d{V) = max \v2{y) - vi{y)\. 
y<^[yo,yi] 

3.1.2 The Conley-Moser conditions 

Let ip : S ^ R he a map and let I = {1, 2, A'"} be an index set for some N £ N. 
Let {Hi}i^f be a set of disjoint m/i-horizontal strips and let {Vi}i(zf be a set of disjoint 
mt,-vertical strips. The Conley-Moser conditions on ip are as follows: 

Condition 3.1.1. ^ mynih < 1- 

Condition 3.1.2. i/; maps Hi homeomorphically onto Vi (i.e. ip{Hi) = Vi) for each 
i G I. Moreover, the horizontal boundaries of Hi map to the horizontal boundaries of 
Vi and the vertical boundaries of Hi map to the vertical boundaries ofVi. 

Condition 3.1.3. Suppose H C Uie/ ^« ^"^ ^''^ m^^-horizontal strip and let 

H^ = ip-\H) n Hi. 

Then Hi is an rrih-horizontal strip for each i £ I and d{Hi) ^ nfid{H) for some < 
rih < 1. 

Similarly suppose V C Uie/ ^ ^"^ '^^ rriy -vertical strip and let 

Then Vi is an m^-vertical strip for each i £ I and d{Vi) ^ nyd{V) for some < n„ < 1. 

The result we will use is the following. 
I 1 

Theorem 3.1.4 flMoser fl973 )). Suppose tp : S ^ R satisfies Conditions [3Al[ [TTE 
and lS.l.M Then ip has an invariant Cantor set A on which it is topologically conjugate 
to a full shift on N symbols, i.e. the following diagram commutes: 
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A " A 



0" 



where </> : A — > is a homeomorphism and a : is i/ie s/ii/t map on the space 

of symbol sequences, defined in section [2.1.SX 

We remark that the restriction of to A has the Bernoulh property. 

3.2 Construction of the strips 

We now construct strips satisfying the Conley-Moser conditions of the previous section. 
Our first task will be to define a certain quadrilateral M C S. It will transpire that the 
images and pre-images of M (with respect to the linked- twist map oF^) have certain 
convenient properties; in particular they contain the horizontal and vertical strips we 
require. 

3.2.1 Construction of the quadrilateral M C S 

Figures [3.21 and [3.31 illustrate the construction of M, which we now describe. Recall our 
notation for the manifold R C T^, established in Section [1.3.H in particular we have 
a 'horizontal' annulus P with boundaries OPq (on which y = yo) and dPi (on which 
y = yi), and a 'vertical' annulus Q with boundaries dQo (on which x = xq) and dQi 
(on which x = xi). 

Consider the portion of the boundary of Q given by dQo D P. This is shown in 
Fig I3.2l fa). Let Si = F~^dQo HP), as shown in I3.2r b). For illustrative purposes we 
have taken j = 2. We observe that Si n consists of j disjoint pieces, each of which 
stretches across S from dQo to dQi. 

One of these pieces has as an end-point the point (xo,yi). Let ffi C Si n denote 
this piece, as shown in I3.2l fc). Similarly we define ff2 C S2 H S, which is derived in a 
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Figure 3.2: Construction of the horizontal boundaries of M. 



similar manner from S2 = F~^{dQi H P) and shown in the same figure; it has as an 
end-point the point (a;i,yo)- 

Analogously, part (a) of Figure 13.31 shows dPo n Q and I3.3l fb) shows its image with 
respect to the map . We have illustrated this using A; = 3. Define Ti = G^{dPQ n Q), 
then Ti n S" consists of k disjoint pieces, each of which stretches across S from bottom 
to top. 

Let fi C Ti n S be that piece which has as an end-point the point [xq^uq). Similarly 
define f2 to be that piece of r2 n which has as an end-point the point (xi,yi); here 

T2 = G^{^Plr^Q). 

Part (c) of Figure [331 shows the quadrilateral M C S, which is bounded by the four 
lines (Ti,(f2,Ti and 7^2- Notice that the boundary consists of two pairs of parallel lines 
and so M is a parallelogram. Finally, we denote by ai C (fi that part of di which is a 
part of the boundary of M, and similarly a2 C (72, ti C ri and r2 C T2- 




Figure 3.3: Construction of the vertical boundaries of M. 
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3.2.2 Intersection of M with its images and pre-images 

In this section we prove results concerning the intersection of M with its image and 
pre- image (with respect to the Hnked- twist map o F^) respectively. The result- 
ing sets resemble but technically are not collections of m^,-vertical and m/j-horizontal 
strips, because they do not stretch completely across S. It turns out that they are the 
intersection of such strips with M. We state the required result as a proposition. 
We recall that the inverse of o is given by F^^ o G^^. 

Proposition 3.2.1. M has the following properties: 

1. M n {G^ o F^{M)) is the intersection of M with {j — l){k — 1) disjoint m^-vertical 
strips. These strips intersect only the boundaries ai and cj2 of M (i.e. they do not 
intersect the boundaries ti or t^). 

2. Similarly, {F^^ oG^^{M))r\M is the intersection of M with {j — \){k — \) disjoint 
nifi-horizontal strips. These strips intersect only the boundaries ti and T2 of M 
(i.e. they do not intersect the boundaries o"i or 02). 

3. The above holds with < m^mh < 1. 
Proof. We prove each part in turn. 



g\f\m)) 



Figure 3.4: Construction of vertical strips crossing S. 

1. Consider F^{M), as illustrated in Figure [3^ a) with j = 2. By the construction of 
M this strip crosses 5 horizontally (j — 1) times and, because the linked-twist map 
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a 




is a homeomorphism, no two of these crossings intersect each other. Part (b) of the 
figure shows o F^{M), ihustrated with k = 3. Observe that {G^ o F^{M)) n S 
consists of (j — 1)(A; + 1) = 4 disjoint pieces. 

The boundaries of M are straight hnes and / and g are affine, so the boundaries 
of these pieces are straight hnes. By definition, the (j — l)(/c — 1) = 2 such pieces 
that cross S completely (i.e. intersect both OPq and dPi) are disjoint mt,-vertical 
strips. 

In Figure I3.6f a) we show the same situation in closer detail and with the original 
set M overlaid. To complete the proof of the first part of the proposition, it suffices 
to show that none of the (j — pieces of {G^oF^ {M))nS intersect (fiUr2) C 
G^ o F^{dP). Assume for a contradiction that this is not the case. Consequently 
we can find some {x,y) G R for which {x,y) G {G^ o F^{M)) Pi [G^ o F^{dP)). 
Then {F~^ o G^^{x, y)) G M n dP = 0, an obvious contradiction. 




Figure 3.5: Construction of horizontal strips crossing S. 



2. The proof of the second part is similar. G^^{M) crosses 5 vertically {k — 1) times 
and these strips are disjoint; see Figure [331 a) . {F^^ o G^^{M)) n S consists of 
+ 1) = 6 disjoint pieces each bounded by straight lines as shown in part (b) 
of the figure. Those pieces which cross S completely are disjoint m/j-horizontal 
strips. 

Now consider Figure [331 fb). In order to show the second part of the lemma it 



52 



suffices to sfiow tliat none of tliese (j + l){k — 1) pieces intersect ai or cj2- This 
follows (as in part 1 of the proof) from the fact that any such point would have 
an image in M n dQ = 0, a clear contradiction. 

3. Last of all we consider the size of m^mh- It is clear from Figure 13.61 part (b) 
that each boundary of each m/j-horizontal strip is an m/i-horizontal curve, and 
moreover that it is a straight line of constant gradient 



yi - yo 



j + Xb- Xa 



with respect to x. Here Xa and Xf, are as shown in the figure. It will be enough 
for present purposes to observe that < Xb — Xa < Xi — xq < 1, thus this gradient 
(which is a suitable value for m/j) satisfies 

yi-yo ^ ^ yi-yo 

. , . <mh< : . 

Similarly from Figure [3T6l f a) we conclude that rriu may be taken to be the gradient 

xi - Xq 



k + yb-ya 



with respect to y. The values i/a and yb are as shown in the figure and satisfy 
< yb — ya < ^- Consequently 



Xi — Xq Xi — Xq 

< ruv < ; 



k + l 
and thus 



(xi - xo){yi - yo) {xi - xo){yi - yo) 

{j+l){k + l) jk 
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□ 




Figure 3.6: Intersection of M with the horizontal and vertical strips. 

Recall our notation that N = {j — 1)(A; — 1) (and our assumption that N ^ 2) and 
/ = {1, 2, N} is an index set. Let {i^jjjg/ be the connected, pairwise-disjoint subsets 
of {F-^ o G-^{M)) n M. Proposition EXU says that for each i G I we have 

Hi = Mn Hi, (3.2.1) 

where Hi is an m/j-horizontal strip. Similarly, let {V^jig/ be the connected, pairwise- 
disjoint subsets of M n {G^ o F^{M)). Then for each i I we have 

Vi = MnVi, (3.2.2) 

where Vi is an m^j-vertical strip. 

3.2.3 Existence of the horseshoe 

In this section we show that the strips we have constructed satisfy the Conley-Moser 
conditions. 
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Proposition 3.2.2. For each i £ I the linked-twist map o maps Hi homeomor- 
phically onto Vi (up to permutation of the Vi). Moreover the horizontal boundaries of 
Hi are mapped to the horizontal boundaries ofVi, and the vertical boundaries of Hi are 
mapped to the vertical boundaries of Vi . 

Proof. o F^ is a homeomorphism of R = P U Q and so it is immediate that any 
one of the A'^ disjoint pieces Hi (a connected component of {F^^ o G^^{M)) n M) 
must map to one and only one of the N disjoint pieces Vi (the connected components 
of M n {G^ o F^{M))), and must do so honi eomorphicahy. Furthermore, elementary 



Armstrong (19831 )) tells us that boundaries must map to 



topology (see, for example, 
boundaries. 

The horizontal boundaries of the Hi are shown in Figure I3.5r b) and their images 
under F^ , shown in part (a) of that figure, are contained within the vertical boundaries of 
G~'=(M). Applying the map G^ in turn, these vertical boundaries become the horizontal 
boundaries of M, which contain the horizontal boundaries of Vi. 

Similarly, the vertical boundaries of the Vi are shown in Figure [3^ bl. The map 
G"'^ takes these into the horizontal boundaries of F^(M), as in part (a) of that figure, 
and these in turn are mapped by F~^ into the vertical boundaries of M. The vertical 
boundaries of the Hi are contained within these vertical boundaries of M and thus the 
result. □ 

At this point we discuss how the Conley-Moser conditions (as given) are not perfectly 
suited to the present task, how we get around this and how else we might have gotten 
around this. 

The Conley-Moser conditions describe quite specifically how horizontal and vertical 
curves and strips are mapped onto each other. We have defined curves as stretching 
completely across S and consequently the boundaries do not map in the required man- 
ner; to remedy this we consider the intersections of such curves with M, as in (13.2. ip 
and (I3X2]1 . 
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Defining curves and strips as we do allows us to determine our estimates on widths 
using the orthogonal {x, y) coordinates on at the expense of then having to intersect 
these curves and strips with M in order that the Conley-Moser conditions are satisfied. 

A different approach would be, as we have suggested, to define the curves and strips 
only on M; the pay-off here is obvious, but we are then forced to adopt new coordinates 
on M and this clearly complicates matters in its own way. We are of the opinion that 
the former is the 'lesser of two evils'. 

Consequently we are forced to adopt the (cumbersome) notation that strips denoted 
with an overbar stretch across 5, whereas strips denoted without an overbar stretch 
only across M . Thus H represents an m/j-horizontal strip whereas H represents the 
corresponding intersection H Ci M. 

Proposition 3.2.3. The proposition has two parts: 

1. Let H he an mh-horizontal strip such that H = H n M is contained in Uiez-^i 
and let 

Hi = {F-^ oG-''{H))nHi. 
Then there exists an ruh-horizontal strip Hi such that Hi = HiH M and 

d{Hi) ^ nhd{H) 

for some < n/j < 1. 

2. Similarly, let V be an my -vertical strip such that V = V D M is contained in 

Vi = {G''oF^{V))nVi. 
Then there exists an -vertical strip Vi such that = H M and 



d{Vi) ^ nyd{V) 
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for some < < 1. 




Figure 3.7: Pre-image of horizontal strips: Part I 

Proof. We prove just the first part, the second being similar. 

Figure [3?71 fa) shows an m/j-horizontal strip H and two strips across M, Hi = H = 
H f] M and H2 (both shaded). Figure ISTTT b) shows the pre-images of M, H = Hi and 
H2 with respect to G'^, illustrated with /c = 3. We observe that the boundaries of each 
are straight lines and that each of G^'^iM) n S, G'^iH) n S and G-^[H2) n S consists 
of {k - 1) m^-vertical curves. Moreover {G-^{H U H2) r^ S) C G^'^iM) n S. 

Recah from Figure[33]^b) that {F^^ o G'''{M)) n S consists of (j + l){k - 1) disjoint 
pieces and that (j — l){k — 1) of these are m/j-horizontal strips. Figure [3^ a) shows 
only these pieces. Notice that o G~^{H) (and similarly F"^ o G^'^{H2)) stretches 
completely across each piece and has straight-line boundaries. In other words, {F~^ o 
G^^(H)) n S consists of (j — l)(k — 1) m/j-horizontal strips. We denote these \ Hi\ , 
SO that Hi = Hif] M. 

It remains to show that d{Hi) ^ nhd{H) for some < < 1. Each Hi is the 
intersection of M with an m/j-horizontal strip of width d{H). Similarly each intersection 
Hi n [F^^ oG^^{Hii)) (for i, i' G /), is the intersection of M with another m/i-horizontal 
strip of width d{Hi). To obtain the result, consider Figure [3^ b): it is clear that 

(i - 1)(A; - l)d{&i) ^ d{H). 
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Simply observe that d{Hi) = d{Hi) and that d{H) = d{H), and one obtains the required 
result with 




□ 



Figure 3.8: Pre-image of horizontal strips: Part II 



We are now in a position to prove the existence of the horseshoe. 

Theorem 3.2.4. Let each of j,k and the product N = {j — l){k — 1) be at least 2. 
Then the toral Unked-twist map o has an invariant Cantor set A on which it is 
topologically conjugate to a full shift on N = (j — l){k — 1) symbols. 

Proof. We show that Propositions 13. 2. H 13.2.21 and 13.2.31 together imply that Condi- 
tions 13.1.11 13.1.21 and 13.1.31 hold. The result then follows from Theorem 13.1.41 

Proposition 13. 2. ll savs that o {M) contains N disjoint m„-vertical strips, which 
do not intersect the boundaries ri or T2 of M. Denote by Hi these strips and by Hi 
their respective intersections with M. Similarly F^^ o G^^{M) contains N disjoint 
m/i-horizontal strips which do not intersect ai or (72. We denote these by Vi and de- 
note by Vi their respective intersections with M. We have < m^mh < 1, satisfying 
Condition EXH 

It is important to notice that HiOVj = HiCi Vj for all i,j £ I, i.e. no intersections 
between horizontal and vertical strips occur outside of M. 
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Proposition 13.2.21 shows that G'^ o acts as dictated by Condition 13.1.21 i.e. the 
Hi are mapped homeomorphically onto the Vi and horizontal (respectively, vertical) 
boundaries are mapped to horizontal (vertical) boundaries. Thus o satisfies 
Condition [3X21 

Finally, Proposition 13.2.31 shows that G^ o F^ acts as specified by Condition 13.1.31 
In particular the pre-image F~^ o G~^{Hi) forms another horizontal strips when it 
intersects with the strips Hi , and each of the new strips has width strictly less than the 
original. Analogous behaviour occurs for the images G'^oF^iVi). Thus Proposition EXl 
shows that G^ o F^ satisfies Condition 13.1.31 □ 



3.3 Uniform hyperbolicity of the horseshoe 

We conclude the chapter with a proof that the restriction of G^ o F^ to the invariant 
set A satisfies the uniform hyperbolicity conditions given in section 12.21 This result 
should not be surprising to us. The toral linked-twist map we consider is seen to be a 
generalisation of the uniformly hyperbolic toral automorphism known as the cat map. 
The non-uniformity of the present system derives from the fact that stretching and 
contraction, the very essence of hyperbolicity, are a conseq uence of points entering and 



returning to S (this feature is highlighted in Wojtkowski's ([19801) proof of the Bernoulli 
property). That we have no lower bounds on this return-time means that any uniform 
hyperbolicity constants we propose can be violated by some trajectory. 

With A the situation is different. By construction the return-time to S is one 
iteration for each point. Let z G A; we have the Jacobian 



D{G'' o F^] 




1 ja 

.kp l+jkal3j 



which is independent of z itself. Denote this Jacobian by A for convenience, then the 
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eigenvalues of A are given by 




It is easily checked that these are real and distinct, and moreover that 



< A_ < 1< A 



-1 



Eigenvectors of A are given by 



v± = 




Define subspaces E^{z) and E'^{z) in the tangent space T^T^ to be the spans of V- 
and respectively. Because these vectors are linearly independent they form a basis 
for the tangent space, i.e. T^T^ = ©i?"(z). The properties required of a uniformly 

hyperbolic system are easily satisfied with C = 1 and A = A_ ; indeed 



and similarly = {G^ o F\z)) . Together these satisfy (f2XT]) . Finally let 

Vs G E'^{z) and Vu G E'"-{z), then 

= WXlvsW = \X":\\\vs\\ and = ||A+"f„|| = |A"|||u,,||, 

where || • || denotes the standard Euclidean norm; these satisfy ()2.2.2p and ()2.2.3p re- 
spectively. 



AE%z) = {Apv- :p£R} = {pA_u_ : p e M} = E' [G'' o F\z) 
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4 The Bernoulli property for a 
planar linked-twist map 



In this chapter we estabhsh the Bernoulh property for the planar Hnked-twist map 
defined in Section I1.3.2I Our starting point is the work of iWoitkowski (198Q 1 who 
proved that such systems can have an ergodic partition. We outline his criteria for this 
in Section 14. 1[ 

A crucial element of our proof is the introduction of new coordinates on the manifold 
A which will enable us to improve Wojtkowski's estimates on the orientation of local 
invariant manifolds. We introduce these in Section 14. 2i They simplify the proof that in 
the case 2 ^ rg < ri ^ ^/7, Wojtkowski's condition for an ergodic partition is satisfied 
and we give this proof in the same section. 

In Section [4.31 we describe the local unstable manifold j'^{w) of a point w = (n, v) £ 
A. Following Wojtkowski's lead it will be convenient to consider the unstable manifolds 
7"(a;), where uj = (r, 9) = Ml^^(w). We define its length and discuss how this definition 
might be extended to Q'^{'y'^ (uj)) . There are some technical considerations but we show 
that for n-a.e. such w our definition does indeed hold. We conclude by showing that 
the length of 0^(7"(w)) grows exponentially with n. 

In Section [4.41 we express the planar linked-twist map G in our new coordinates. We 
then show that a certain tangent cone is preserved by the differential DQ. Finally, in 
Section [53] we are able to give the estimate we have mentioned on the orientation of the 
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unstable manifolds and, following this, a largely geometrical proof that the Bernoulli 
property is satisfied. 



4.1 Wojtkowski's results 



In this section we describe Wojtkowski's (| 19801 ) criteria for the planar linked-twist map 
to have an ergodic partition. 

Let w = {u,v) £ T, and denote by a{w) G (0, vr) the angle at which the segment 
connecting w to (—1,0) meets the segment connecting w to (1,0). Let 



cot a(w) / . - , X 

r? = sup ^-l-l 



where r(w) denotes the Euclidean distance from w to (—1,0). We also denote by c 
the infimum of the derivative of the twist function, which in this case is just given by 
27r/(ri — ro). Wojtkowski proved the following: 



Theorem 4.1.1 (jWoitkowski figsd ll. // 

c>2r] (4.1.2) 

then the linked-twist map Q : A ^ A is the union of (at most) countably many ergodic 
components. 

Wojtkowski's conjecture, made in the same paper, is that under the assumptions 
of Theorem 14.1.11 then : A ^ A has the X-property. This would, by the work of 



Chernov and Haskell (19961 ) . imply that it has the Bernoulli property. 



We discuss the proof of Theorem 14.1.11 briefly. It is easily argued that ;U-a.e. w £ A 
lands in S under iteration of and, furthermore, returns to S infinitely many timesll] 



^One simply notices that those points not satisfying this condition must be rigid rotations around 
one of the annuh, and must have rational angle of rotation, else their orbit would be dense and hit E. 
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Thus for a full- measure set of points we may talk of the return map to T,, o r just the 



return map as we shall usually abbreviate it. Following IWojtkowski (19801 ) we shall 
actually define the return map on MJ^^{Y,) C L rather than on S itself, as follows: 

Definition (First-return map to S). Let {r,9) G M^^{T,) C L. The first-return map 
to S is the map Q^, ■ M^^{Yi) — > M^^{Y,) given by 



92 = o o M+, 



where i is the smallest (strictly) positive integer for which Q^{M^{r,9)) G S. 

For {r,e) G M+^(S) let /3i = dr, P2 = d9 give coordinates in the tangent space 
T(f, g)L and define the cone 



U{r,e) = \{P,,P2):^^>^ 



Wojtkowski establishes that U is invariant under, and expanded by, the derivative 
We illustrate the situation in Figure I4.1[ More precisely, define the cone field 

U+= [j Uir,9) 
(r,e)eM-i{s) 



and let || • || be the norm in T(^r,e)L induced by the Riemannian metric, i.e. ||(/3i,/32) 



YT^i+r^^ . We have the following: 



Proposition 4.1.2 (IWojtkowski (1980 )). Z)0e({7+) C {/+. Furthermore there is a 



constant A > 1, independent of {r,6) or (5, and for vectors (3 G we have \\DQy.P\\ ^ 



For a proof see the original paper or 



Sturman et al. (20061 ) 



From the nature of the twist function (in particular the strict monotonicity) one easily infers that such 
points are contained within a set of measure zero. 
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dr 



dr 



d^ 



e 




Figure 4.1: The invariant expansive cone U C T^^L is shown in the left-hand figure. In 
the right-hand figure is the image of the cone under the differential map DQy. (dark- 
shaded) with the original cone (light-shaded) included for comparison. Observe how 
the cone is mapped into itself and vectors within it are expanded. 

4.2 New coordinates for the manifold A 

At the heart of our proof of the Bernoulli property for is a new coordinate system. In 
this section we introduce these coordinates and with them prove that if the condition 



is satisfied, then the linked-twist map has an ergodic partition. Theorem 14.1.11 savs 
that this amounts to proving that the condition (j4.1.2p is satisfied. 

We use the notation = [— vr, vr] with opposite ends identified, TZ = [ro,ri] and 



4.2.1 Construction of the new coordinates on 

We introduce new coordinates {x{u,v),y{u,v)) for the manifold A. The key feature is 
that if {u, v) G S+ then x{u, v) is given by the Euclidean distance from (u, v) to the 
centre of annulus Aj^ and y{u,v) is given by the distance from {u,v) to the centre of 
annulus ^-1^ Thus for {u,v) G S+ we have 

^Coordinates such as these are commonly caUed two-centre hi-polar coordinates, however we do not 
adopt this name because we will not extend the coordinates to all of A in this manner. 



2 ^ ro < ri ^ \/7 



(4.2.1) 



n = [-ri,-ro]. 
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(x,y) = (^y/{l + u)^ + v\ ^(l_^)2 + „2^ 



(4.2.2) 



Similarly on S_ the magnitudes of x and y are determined in this way, but one of the 
coordinates assumes a negative value so that points are uniquely defined. Unfortunately 
it is not useful to extend x, y to all of A in the same manner so we will take an alternative 
approach. We begin by defining the coordinates on A^. 

Let (u, v) £ A-^ and let (r, 6) = M^^{u, v) L = TZ xS^. Setting x = r satisfies the 
first part of (14.2.20 because r is the Euclidean distance from {u,v) to (—1,0). However 
the coordinate will not in general give the Euclidean distance to (1,0). 

We now wish to define a homeomorphism ^' : 7^ x ^ 7^ x so that 



for some function ip : TZ x S which is to be determined. It suffices to define a 

homeomorphism -0 : 7^ x [0, vr] — > [0, vr], for which '(/'(r, 0) = 0, '(/'(r, vr) = vr and with the 
extra condition that 



so that ()4.2.2p holds. Once we have done this, we will extend our coordinates to A-.. 
Figure 14.21 should help the reader to keep track of our construction. 

Let {u,v) G S+, which is one of the light-shaded regions in Figure [121 We require 
that ip{r, 9) = y^(l — u)'^ + v'^. Expressing u, v in terms of r, 6 gives 



{x,y) = ^ o M^^{u,v) 



is as required. It is clear that ^ will need to have the form 



i>{r,-9) 



(4.2.3) 




(4.2.4) 
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V By 



Figure 4.2: The region yl+ C M^, illustrated in the three coordinate systems. Left- 
to-right: Cartesians (u^v) in the plane; polars (r, ^) G x and new coordinates 
{x^y) G X S-*^. Shading indicates the three regions for which ^ takes different forms, 
as explained in the text. 
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This defines V (and thus the new coordinates {x,y)) on the domain S = S+ U S_. 
It is clear from the construction (of two-centre bi-polar coordinates) that this is a 
homeomorphism. 

Before we continue it wih be convenient to determine the inverse to ()4.2.4p onTZxTZ 
, i.e. the function ^"""^ : TZ x TZ ^ Af^^($]) such that 

'ilj{x,'4!~^{x,y)) = y. 

We claim that the function cos^^ (r(x,y)), where T{x,y) = ^ (x^ — + 4), has this 
property. Indeed: 

■0 (x, cos""^ (t(x, y))) = sgn(y ) x'^ — 4xr (x, y) + 4 

= sgn(y) Y^x^ — + — 4 + 4 = y. 

We now define new coordinates on j4_|_\S, which consists of two disjoint pieces (recall 
that we are considering only the 'top half as illustrated; we will deal with the bottom 
half using the symmetry given by (|4.2.3p ). We deal with each separately. Our approach 
to this is very simple though perhaps this simplicity can be lost amongst the equations 
we provide. To counter this let us describe heuristically what it is that we are doing 
first. 

Consider Figure 14.21 again and look at the middle figure, in (r, 6) coordinates. We 
have constructed a homeomorphism from M^^{T,j^), which is the upper, light-shaded 
area, to the corresponding light-shaded area in the right-hand figure. Now we want to 
do the same for the dark-shaded area immediately below M_^^(S+), that is to create a 
homeomorphism from it to the corresponding dark-shaded area in the right-hand figure. 
The simplest way to do this is to leave the points such that 9 = invariant and 'stretch' 
the others, 'upwards' as illustrated, whilst leaving the r-coordinate unchanged. 

More formally: let {u,v) be a point in that part of which lies inside the 
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annulus A_. This is the dark-shaded region in Figure HT2l M_^_^{u,v) has coordinates 
{r,6) such that r £ TZ and 6 £ [0, cos~^(T(r, ro))]. For each r £ TZ we define on 
{r} X [O,cos"i(r(r,ro))] by 

^(r, 6*) = ro6'/ cos"^(r(r, tq)). (4.2.5) 

Finahy let (n, denote a point in that part of which hes outside of A^. 

This is the unshaded region in Figure Here AI^^{u,v) is given by coordinates 

(r, 0) £ TZx [cos~-'^(r(r, ri)), vr]. For each r £ TZ we define ^jJ on {r} x [cos^^(r(r, ri)), vr] 

by 

^(r,^)=n+ ^-^""/:f'^^^,^. (vr-n). (4.2.6) 
TT — cos ^(T(r, rijj 

This, together with the symmetry requirement (I4.2.3P completes our definition of x{u, v),y{ 
on A^. To summarise: 

Definition (Coordinates {x,y) on A^). Let {u,v) £ A^, then {x,y) = ^' o AI^^{u,v), 
where ^(r,6) = (r,ilj{r,6)) and tp takes one of the forms 1^4. 2. 4% li4-2.5 ) or Ii4.2.6^ as 
described. 

4.2.2 Construction of the new coordinates on A_ 

There is a slight problem to overcome in extending the coordinates to A-. Because 
of the geometry of A, the most natural approach (i.e. straight-forward symmetry) will 
lead to us having two different expressions for the new coordinates on S_ . Our strategy 
is to pursue this naive approach anyway and deal with the problem when it arises; in 
acknowledgment of this we call the coordinates on A- 'temporary' for the moment. We 
hope that Figure 14.31 will help the reader. 

We will define temporary new coordinates {x',y') on A^ before extending {x,y) to 
that domain. We introduce the functions l, : M? ^ given by 
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i^ix^y) = {y,-x) and l ^{x,y) = {-y,x). 

With ^ as already defined, on A- we define {x',y') = l^^ o ^1/ o AIZ^ {u,v). Notice that 
MZ^ expresses A- in polars (albeit centred at (1,0) and with polar angles shifted by 
vr), so that the restriction of ^' to MZ^{A_) is a homeomorphism, just as the restriction 
of ^ to M^^(^+) was. For {u,v) G S-t we have 

{x',y') = L-'^ o^i o Mz'^{u,v) 

= i-^ (V(i-^^)' + ^^^tV(i + ^)' + ^') 

Comparing the coordinates {x,y) and {x',y') where their domains overlap we see that 
{x\y') = ±(x,y) on T,±. See Figure We thus extend the coordinates (x,y) to A^ 
as follows 

Definition (Coordinates {x,y) on A^). For {u,v) £ yl_\S_ define new coordinates 
{x, y) = (x', y') = o MZ^{u, v). For {u, v) G S_ define {x, y) = —{x',y'). 




Figure 4.3: The manifold A, illustrated in its native Cartesian coordinates and in the 
new coordinates {x,y) G x S-*^. Notice that there are two distinct representations 
of S_ deriving from the coordinate transformations. In (x, y)-coordinates this image 
is given by 7^ x —TZ (toward the bottom-right of the right-hand figure), whereas in 
(x', y')-coordinates it is given by —TZ x TZ (toward the top-left of the right-hand figure). 
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4.2.3 The ergodic partition 

The new coordinates simplify our proof that the planar linked-twist map Q : A ^ A 
defined with radii 2 ^ tq < ri ^ -v/7 has an ergodic partition, as defined within 
Theorem 12.2.51 This amounts to the following: 



Lemma 4.2.1. Let 2 ^ tq < ri ^ ^/7. Then condition \4^.1.2 ) is satisfied. 



Proof. By symmetry it is enough to show that the condition holds on Moreover 
the condition is implied by 

sup cotQ!(?i;) < — . (4.2.7) 
i«ei;+ v7 — 2 

Let {x,y) = o AI^^{w) give w in the new coordinates. The angle a appears in a 
triangle in which its adjacent sides have lengths x and y, and the opposite side has 
length 2. The law of cosines says that 

cosa = ^^i^-^. (4.2.8) 
2xy 

The partial derivative of (|4.2.8p with respect to x is given by 

dx 2y 2x^y ' 

and, using x,y £ [2,\/7], we calculate that l/2y E [l/2\/7, 1/4] and (y^ - 4)/2x^y e 
[0, 3/16]. It is easily checked that 1/2^/7 > 3/16 and so the derivative is always positive. 
Consequently cos a is an increasing function of x and thus a is a decreasing function of 
X. By symmetry a is also a decreasing function of y. 

Combining these facts with ()4.2.8p we find that a £ [cos^"'^(5/7), 7r/3]. Recall that 
cot is positive and decreasing on (0,7r/2), then 
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„i 5 5^6 

sup cot a = cot mi a = cot cos - = , 
weT,+ »«ge+ 7 12 



and (I4.2.7P is seen to be satisfied. 



□ 



4.3 Growth of local invariant manifolds 

In this section we describe the nature of local unstable manifolds and their images 
under Q, and define what we mean by the ir length. We show that this length diverges 



on iteration of the map. We remark that 
are true but did not give a proof. 



Woitkowski (1980 ) stated that these results 



4.3.1 Nature of local unstable manifolds 

At /i-a.e. point (u, v) € T j there is a positive Lyapunov exponent. This follows easily 
from Wojtkowski's (|l980|) construction of the in variant expansive co ne U, although he 



does not prove it explicitly (his work pre-dated iKatok et al. (19861 ) so he was unable 



to use their extension of Pesin theory, a,nd thus had no reas on to discuss Lyapunov 



exponents). The details r nay be found in 



Theorem [223] (due to 



Sturman et al. (20061 ). however 



Katok et al. (19861 )) tells us that as a consequence, associated 



to each such point is a local unstable manifold. Let w = {r,9) = M_^^(u,v). The local 
unstable manifold is denoted 7"(w) and has the form 



(4.3.1) 



for some e > 0. Here, i?"(0, e) is the open e-neighbourhood of the origin in the unstable 
subspace E'^{oj) C T^^L and 0" : i?"(0,e) — > E^{uj) C T^^L is a smooth map satisfying 
0"(O) = and D</>"(0) = 0. 

We prove that the planar linked-twist map Q : A ^ A has the Bernoulli property by 
demonstrating that the 'strong' form of the condition (j2.2.13p is satisfied. We recall that 
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/x-a.e. point returns infinitely many times to S and that we are defining local invariant 
manifolds in M_jr^($]) rather than in S itself. The condition is that for /x-a.e. Wjw' gT, 
and for all sufficiently large integers m and n (depending on w,w') we have 

e" (M+ (7"(M+iH))) n Q-^ (M+ (7^(M+i(u;')))) / 0. (4.3.2) 

There are two key components to our proof that the condition (|4.3.2p is satisfied. 
The first, which we deal with in this section, is that local unstable manifolds grow under 
iteration of 0. We make this notion precise once we have defined what exactly we mean 
by their length. The second component is a useful characterisation of their direction 
or orientation. The latter is our principle motivation for having introduced the new 
coordinates and we return to it in the following section. 

4.3.2 Length of local unstable manifolds 



We recall some terminology from differential geometry, taken again from lDo Carmo (19761 ) 
if M is a smooth manifold and / C M an open interval, then a smooth curve is a smooth 
map a : I ^ M. The image 7 = a(/) C M is referred to as the trace of the curve. 

Definition (Length of a smooth curve; length of its trace). Let M be a smooth manifold, 

I C R an open interval and a : I ^ M a smooth curve. We define the length of the 
curve a to be 

length{a) = J \\Dat\\dt, (4.3.3) 

where Dat denotes the derivative (i.e. the Jacobian) of a evaluated at t €z I and \\ ■ \\ 
denotes the Euclidean norm in T^i^i^M . We define the length of the trace 7 = a{I) to 
coincide with the length of the curve a. 

We remark that it is perhaps more common in the differential geometry literature 
to see the expression da / dt in place of Dat but it will be more convenient for us to use 
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the latter expression. We also comment that this definiti on of length can b e shown to 



agree with our geometrical intuition. For more details see 



Do Carmo f 19761 ). 



Now let to = {r,9) G M^^{T,) C L be such that 7"(a;) exists and is of the form 
(j4.3.ip . The unstable subspace and thus the unstable manifold are one-dimensional, so 
there is an open interval / C M and a diffeomorphism a : / — > M^^{A) (i.e. the map 
is diffeomorphic between I and its image in M~^(A)) such that a{I) = 7"(w). The 
diffeomorphism a is a smooth curve of which ^^{uj) is the trace, so the length of 'J^{uj) 
is defined as above. 

Now let n G N and consider the image of 7"(w) with respect to G^- In general 
0^ oa : I ^ M^^{A) will not be a smooth curve so it is not immediately clear that the 
length of (lu)) is defined. Suppose however that we are able to determine that 

o a is piecewise smooth, in the sense that I decomposes into a countable (at most) 
union of intervals 

oo 

I = {io,k)'j\J[ih,ih+i) (4.3.4) 

h=l 

and that the restriction of 0^ o a to any interval (i/i,i/i+i) is smooth. In this case our 
definition of length applies to each restriction and we naturally determine the length of 
o a by summing over them. We illustrate the situation in Figure 14.41 
Determining that such a partition of / exists for a 'typical' point oj G M^^(I]) is the 
main technical difficulty in proving that the length of its local unstable manifold grows 
under iteration with Q. Our ambition for the remainder of this section is to prove the 
following theorem: 

Theorem 4.3.1. Let n G N and let X > 1 be the constant given by Proposition 

For fx-a.e. u; G S, letu) = M^'^iw), then the lengths ofj^{uj) and @y^{^'^{uj)) are defined 

and 

length {Qli-i'^iuj))) ^ X"" length {j"" (uj)) . (4.3.5) 
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I 







eg(7"(^)) 




eg 



a{I) = 7"(a;) 



eg°«(/i) 



Figure 4.4: The length of local unstable manifold 7"(w) is well-defined because there 
is a diffeomorphism a of an open interval / and a{I) = 7"(w). If / decomposes into a 
countable (at most) family of intervals on which 0^ o a is differentiable then the length 
of @^{'y^{u:)) is well-defined also. 

4.3.3 Proof of Theorem |4XT] 

Our discussion of the length of piecewise smooth curves motivates us to formulate the 
following proposition: 

Proposition 4.3.2. For any n E N, for fi-a.e. G S and for uj = Ml^^{w), the set of 
oj' G 7"('j^) at which DQ'^ does not exist is at most countable. 

In proving Proposition 14.3.21 we will use two lemmas. We state and prove these 
lemmas first, then prove the proposition and last of all prove the theorem. We recom- 
mend that the reader skip forward to the proof of the theorem and refer back to the 
proposition and lemmas in that order. 

Our first lemma is easily stated and proven but we do not know of it anywhere in the 
literature. Loosely speaking, it says that when the traces of two continuous, injective 
curves intersect there are only the two possibilities shown in Figure [^31 

Lemma 4.3.3. Let I, J C M be connected open intervals, let M be a metric space and 
let a : I ^ M and (3 : J ^ M each be continuous curves of finite length. Either there 
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a{I')=P{J') «(/) 



(3{J) 



Figure 4.5: Illustration of Lemma [4.3.3t when the traces of continuous, injective curves 
a : I ^ M and (3 : J ^ M intersect, the set of intersections is either at most countable, 
as in the left-hand figure, or it contains the images of open intervals I' C I and J' C J 
such that «(/') = P{J'), as shown in the right-hand figure. 

are open intervals I' C I and J' C J for which a{I') = (3{J'), or the set a{I) Pi (3{J) is 
at most countable. 

Proof. Assume that there are no intervals /' and J' having the property stated, and let 
to, ti £ I he any two distinct points such that a(to) H /3( J) 7^ and a(ti) n (3{J) / 0. (If 
there do not exist such points to,ti then, of course, we are done.) By the 'no intervals' 
assumption, there must be some t G {to,ti) for which a{t) n /3(J) = 0. Set 



where d denotes the metric on M, and where the distance from the point a{t) to the set 
/3(J) is defined, in the usual way, as the infimum over distances d{a{t),p) for p £ f3{J). 

Now, /3( J)) is a continuous function of s G /, so there is some 6 > and an 

open neighbourhood {t - 5,t + 6) C (to, h) C I such that {t - 6,t + 6) D (3{J) = 0. By 
this argument, the number of intersections Oi{I) H /?( J) can exceed the number of such 
open neighbourhoods by at most one and there can be at most countably many disjoint 
open intervals in the interval I. □ 

Our second lemma concerns the nature of the set of points at which DQ^ is not 
differentiable. 



d{a{t),/3{J)) = e > 0, 
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Lemma 4.3.4. Let n S N. The set of points in M_^_ {A) at which is not differ- 

entiable is given 6y M_^^(si U {(-1,0)}), where 

i-l 

s^=\Je-^{dA+U<^-\dA^)) , (4.3.6) 

h=0 

for some i E N. Moreover M^^(sj) is the union of a single point with a finite union 
/3/i( J/i), where N £ N and for each h = 1,2,..., N the set Jh C M. is an open 
interval and j3h '■ Jh ^ Ml^^{A) is a continuous, injective curve of finite length. 

In Figure 14.61 we illustrate the set sq of non-difFerentiable points. 




Figure 4.6: The manifold A showing those points at which Q is non-differentiable. Red 
lines indicate the set so of non-differentiable points, black lines indicate other boundary 
points. 

Proof. We prove the statement in the case n = 1, the proof of the general case is no 
more difficult. RecaU that DQj^ : M+^(S) M^^(S) is given by o 6* o M+ for 
some positive integer i. We investigate whether the derivative 

D(es). = Z)(Af-i)e.(M+H)^0V+M^(^+)- (4.3.7) 
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exists for a given uj E {A). For an affirmative answer we require that each of 
the three derivatives on the right-hand side of (j4.3.7p exist. First note that M+ is 
differ entiable on A. 

Q = r o <I) where <I> is differentiable except on dAj^ and F is differentiable except 
on dA_. Thus we require that M+(u;) ^ dAj^ and that ^{Mj^{uj)) ^ i.e. that 

M+(tj) ^ So = {dA+ U ^'^{dA^)). 

is differentiable if M+(u;) ^ sq and e(M+(cj)) ^ sq, i-e- if M^{uj) ^ soUe-^(so). 
Continuing inductively we see that 0* is differentiable except at those Mj^{u}) G U/iZ=o 

Lastly we observe that M^"^ : A x S-*^ is differentiable except at the point 

(—1,0). Combining these conditions we conclude that the set of points in M^'^{A) at 
which Z)02 is not differentiable is given by M+^(siU{(-l,0)}) for some i G N, with Si 
as defined by ()4.3.6p . 

To see that the second statement holds, observe that (a pair of circles) is the 
trace of two smooth curves. Applying <I>^^ will leave parts of these circles invariant and 
skew those parts which cross yl+ around that annulus. See Figure 14.61 This process 
stretches each circle (by a finite amount) and means that the curve whose trace is 
(^~^{Aj^) is not differentiable at the points which map to dAj^. There are only finitely 
many such points however, so this curve is a finite union of smooth pieces. We proceed 
inductively, observing that each of the maps <I>^^,B^'^ and M_^^ will introduce only a 
finite number of non-differentiable points at each step, to obtain the result. □ 

We are now in a position to prove the proposition. 

Proof of Proposition \4.3'^ Let z S N and let Xi C M^^(S) contain precisely those 
uj for which 7"(cj) H M^^{si) is uncountable, where Si is defined by (j4.3.6l) . The set 
M^^{si) consists of those points at which DQj^ is non-differentiable, so Xi consists of 
those points uj G Af^^(S) for which we might not be able to determine the length of 
Gg(7(a;)). We will show that /u(M+(Xi)) = 0. 
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Lemma 14.3.41 savs that if 7"(w) exists for some uj E M^"'^(S), then the set 7"(w) H 
M^^{si) is a finite union 



N 

[ja{I)r\(ih{Jh). (4.3.8) 
h=i 

where / and each are connected open intervals of M and a : I ^ M^^{A) and each 
f3h '■ Jh ^ M^'^{A) are continuous curves of finite length. Now assume that the set 
(j4.3.8p is uncountable, i.e. lo E Xj. Then there is at least one h E 1,2,...,A^ for which 
a[I) n (3h{Jh) is uncountable. Fix such an h. 

Lemma 14.3.31 says that there are intervals I' C I and C Jh and that a{I') = 
Ph{J'h)- This is significant because it says that 7"(u;) coincides with M^'^{si) on a 
subset of positive length. The total length of Mjp^(sj) is well-defined and finite, so if 
{'^l^leL is any collection of points uJi E Xi such that ^"^{uji^) n ^^{uji^) = for li / /2, 
then L is at most a countable set. 

Our strategy is to find such a countable set that covers X^. To that end we recall 
the familiar global unstable manifold of u; E given by 



oo 

1=1 



(See for example 



Barreira and Pesin (2002 ).) It is easily shown that uj E W'^{ijJ ) if and 



only if uj' E ^{lo) and that otherwise we have ^{lj) n W^{uj') = 0. Thus for uj ^ Xi 
there is a well-defined equivalence class 

[uj] = [uj' £Xi:uj' ^Wiuj)] . 

We can clearly cover Xi for some set of distinct elements {[i^i]}ieL and by our previous 
reasoning any such L is at most countable. This is sufficient to give the result: 
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oo 



//(M+(X,)) ^^i[M+[ (j[u;i] ^ ^/x(M+(M)) = 0. 



\l=i / / 1=1 



□ 

For any given n G N we can use Proposition 14.3.21 and our definition of tlie lengtli 
of the trace of a smootli curve to define the length of 0^(7"('-^)) for a set of uj whose 
M+-image has full /x- measure in S. This enables us to show that this length grows 
exponentially with n, as in the following theorem. 

Proof of Theorem \4- 3. 1\ Proposition 14.3.2] says that for /x-a.e. if G S and for lo = 
Mj^^{w) we have Q^{'~f^{uj)) = o oi{I) where / C M is an open interval and 
a : I ^ AI^^{A) is a smooth curve. Moreover there is a countable (at most) parti- 
tion of / as given by (j4.3.4p and the restriction of 0^ o a to each open interval {ih, ih+i) 
is a smooth curve. The length of each smooth curve, by our earlier definition, is given 

by 

length(eS o a|(,,,,,^,)) = / \\DeiDat\\ dt, 



and so the length of 0v o a is given by 



length(e£ o a) = ^to J-^" \\D@lDat\\ dt 



= X^J,\\Dat\\dt 
= A"length(a), 

where N £ {oo}. Here the second line holds because of Proposition 14.1.2] and the 
fact that Dat £ U C T^^f-^L for each t £ I. □ 



79 



4.4 A new invariant cone for G 



In this section we express the map @ : A ^ A in the new coordinates developed in 
Section 14. 2[ Once we have done this we will introduce a new cone in the tangent 
space, which is invariant under the map. The purpose of this construction is to give an 
improved estimate on the orientation of local unstable manifolds. This, combined with 
the growth already established will enable us to show that condition ()4.3.2p is satisfied. 
We will do that in the following, final section. 

4.4.1 The map O expressed in the new coordinates 

We will use notation reminiscent of that we have used for other linked-twist maps. 
We begin by giving some notation for the manifold A when transformed into the new 
coordinates. Let 

R= {^^^o M-^{A+)} U {i"^ o o Mr^(A-\S_)} . 

The set C is shown in Figure HTTl f a). There is a one-to-one correspondence between 
points in R and points in A. Let F : R ^ R denote the map ^ : A-^- A^ in the new 
coordinates. Recall that $ = o A o M^^ , so 

F =^o o $ o M+ o ^r-l 

= ^ o A o VI/-1. 

We observe that F is a homeomorphism of R. 
Now let 

i?' = o M_^^(y4+\S_)} U {r^ o^o Mr^(A_)} . 

This is illustrated in Figure HTTl fb). Again, there is a one-to-one correspondence between 
points in R' and points in A. We let G : R' ^ R' denote the map F : A- — > A^ in the 
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Figure 4.7: The manifolds i?, i?' C = §^ x in parts (a) and (b) respectively . 
Each is in one-to-one correspondence with A, the difference between the two being that 
E_ is represented differently in each. The manifolds have the property that i*" is a 
homeomorphism of R and G is a homeomorphism of R' . 

new coordinates, giving 

G = o V]/ o o O i 
= o F-^ o i. 



We observe that G is a homeomorphism of R' . 

Now consider the linked- twist map 6 expressed in the new coordinates, that is the 
composition H = GoF, which we wish to represent as a map of R into itself. H will be 
a homeomorphism as it is the composition of two homeomorphisms. Recall that when 
we defined the map G = F o $ we changed coordinates between the two mappings; here 
we do likewise. The difference is that this is now only necessary when the trajectory 
lands in (the new representation of) S_. Moreover the coordinate transformation is 
given simply by {x, y) i— > {—x, —y). 

Let us make this algorithm more explicit; recall our notation that TZ = [ro,ri] 
and —TZ = [— ri,— ro]. Let {x,y) € R, and first apply the map F. If and only if 
F{x,y) e TZ X —TZ, reflect both coordinates through the origin. We are now in R'. 
Apply G, i.e. apply t, followed by F~'^, followed by Finally, if and only if we are 
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now in —TZ x 7^, reflect coordinates through the origin again. 

Rather than include expressions both for reflections through the origin and the 
functions we are able to combine the two. We have H : R ^ R given by 



H = n-^ oF-^ o^loF, 



where the function Q is given by 



{x,y) = < 



if {x,y) enx^T^ 



L otherwise. 



Let 5 C -R be the image of the 'intersection region' S, i.e. S = {TZ x Tt) U {TZ x —71). 

For z G 5 we define the return map Hs : S ^ S, defined completely analogously to the 
return map Os- 

We now study the derivative of H. Let -Di, D2 denote the usual differential operators. 
We have 



= *o {x,ip~'^{x,y)) 

= ^' (x, tlj^^{x, y) ± c{x — ro)) 

= (x, ^ [x, ^"^ (x, y) ± c{x -ro))) . 



To simplify the expression we define 



f±ix,y) =ip{x,'tp ^{x,y) ±c{x - ro)) 



then the Jacobians of F^^ are given by 
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We introduce further notation y± = y±{x, y) = ip ^{x, y) ± c{x — rg), then 



Dif±{x,y) = Di'4){x,y±) + D2^{x,y±)[Di'4) ^{x,y)±c] 
D2f±{x,y) = D2'4>{x,y±)D2'4>~^{x,y). 



(4.4.1) 
(4.4.2) 



We use these derivatives to prove a result for DH. Let 61 = dx, 62 = dy give coordinates 
in the tangent space T^T^ to a point z = (x, y) G R, and define the cones 

C{z) = {(61, 62) : 6162 > 0} , C{z) = {(61, 62) : 6162 < 0} 

The cone C is illustrated in Figure [48l We define the cone fields 



C+ = U C{z\ = U C{z). 



dy 



-dx 




Figure 4.8: The invariant cone C C T^R is shown in the left-hand figure. In the right- 
hand figure is the image of the cone under the differential map DHg- The fact that C 
is invariant under this differential is immediately implied by Proposition 14.4.11 Notice 
that we do not claim that this cone is expanded by DHs- 

Proposition 4.4.1. Let tq = 2 and ri = \fl . Then DF and D{Q^^ o F^^ oil) preserve 
the cone field C+. 



The remainder of this section is devoted to proving Proposition 14.4.11 We remark 
that we have had to give explicit sizes for the annuli and so our result is not as general 
as one might hope for. We discuss this more in Chapter [6l 
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4.4.2 Proof of Proposition 14.4. 1] 

Both DO, and DQ^^ map the cone C into the cone C and vice versa. Consequently it 
suffices to show that DF preserves (7+ and DF^^ preserves C_. Let z = (x, y) E A, let 
{bi,b2) £ TjT^ and define 



bx 

biDif±{x,y) +b2D2f±{x,y) 





1 = DF^' 1 


fb,\ 


H 






\b2) 





We have 

^ = ^i/±(x,y) + ^^2/±(a:,y). 
b\ bi 

So it is enough to show that for every (x, y) G 7^ x §^ we have 

±Dif±{x,y)>0 and D2f±{x,y) > 0. (4.4.3) 

In light of ()4.4.ip and (|4.4.2p . the condition (I4.4.3P will follow from suitable bounds on 
the derivatives of ■0 and its inverse. Determining such bounds is our strategy for this 
proof; unfortunately, some extensive calculations will be unavoidable. 
Recall the real- valued function r, defined on ^' o M7^($]) = 7^ x 7^ by 



At 



It will play a significant role. Partial derivatives of r are given by 



^) = i^P + ^'- 4) and ^,r{r,e) = -l, 



so T is an increasing function of r and a decreasing function of 9. This observation is 
enough to prove some bounds on certain functions of r, which we will use in the rest of 
this section. We collect these in the following lemma, though we omit a formal proof. 
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Lemma 4.4.2 (Properties of the function r). Let r denote T(r, tq). Then 



T G 



1 V7 
2' 4 



1 r 1 ^ 
' 2-r=4' 



3 VS 
4'^ 



cos 



-1 



5 TT 



Alternatively, let t denote r(r, ri). T/ien 



r G 



8' 7 



2 3 
3' 4 



s/TT7J.(li,i) 



cos r G 



7 3 
6' 2 



It will be convenient to denote by that part of which lies inside the annulus 
j4„ and by A'^ that part of A^ lying outside of A_. Moreover it will be necessary to 
determine the values 9 may take on each part. On for a given r G TZ, this range is 
given by [O,cos~^ r], where r denotes T(r, ro). Using Lemma 14.4.21 we calculate 



1 vr 

sup 9 = supcos"^ T = cos~^ inf r = cos^^ o ~ "q ' 
(r,e)eA' ^•e7^ ^-e??- 2 3 



Similarly, for a given r £ TZ the range of 9 such that (r, 0) G is given by [cos ^ r, vr], 
where r = r(r, ri). We have 



inf = inf cos ^ r = cos ^ sup r = cos ^ > - . 
" - " ^ 7 6 



(r,6»)eA° re7^ 



We now prove some bounds on the derivatives of ip. 



Lemma 4.4.3. Diip G [O, |) . 



Proof. We consider the three cases separately. Let {r,9) G A\, then 



ro0 



'1 _ 

.2 ry 



(cos ^ r)^ Vl — ' 



where r = r(r, tq). The numerator equates to 9/2, so has range [0, vr/G]. Using 
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Lemma 14.4.21 the denominator lies in [25/54, 7r^\/3/18]. Thus on A^:^ we have Diip £ 
[0,97r/25] C [0,7/6). 

Next let {r,0) G yl° , then 



(vr — cos ^ t)'^\/1 — t"^ 



where r = r(r, ri). A lower bound of is attained when 9 = tt; using Lemma 14.4.21 we 
see calculate that 3/10 is an upper bound, so here Ditp G [0,3/10). 
Lastly let (r, 0) G S, so that 



„ / r — 2 cos ( 



where ^|J = 'il){r,9) is as in ()4.2.4p . On recall that cos9 = T[r,ij))^ so the numerator 
is given by 

^^ + 4) = --; + 



- (H - + 4) = h — . 

2r ^ ^ ^ 2 r 2r 



r — — \r 



By construction, ip{r,e) E [r^n] for (r,6') G S+. Also ^ r/2 - 2/r ^ 3^7/14, ah of 
which means that Dii) G [^/7/7, 5^/7/14]. □ 

Lemma 4.4.4. D21P e 

Proof. As before there are three cases. If (r, 9) G then 

D2ip{r, 9) = ro/ cos^"*^ r. 

where r = r(r, ro). Using Lemma [4.4.21 we easily deduce that D21IJ G [4/3, 12/7]. 
If (r, 9) G then 

D2ip{r,9) = (vr — ri)/(7r — cos~^ r). 
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where r = r(r, ri). Lemma 14.4.21 allows us to deduce the bounds -D2V' ^ (1/^) 1/3)- 
Finally if (r, ^) G S+ then 



A lower bound for sin 9 is calculated as follows: 



inf sinO = inf sin (cos~^ T(r,rQ)) = inf \/l — [rfr, rn)!^ = -. 
(r,6l)eE+ reTe ^ ^ rsTe 4 



Trivially, we also have upper bound 1. From here we easily obtain the bounds D21P € 
(3^7/7,^7). □ 

It is easily verified that the inverse to '0 on S (where -0 is given by (I4.2.4j) ) is given 

by 

ip-\x, y) = cos-^ (t(x, y)) . {AAA) 
The inverse to on 7^ x [0, ro] (the inverse to (|4.2.5p ) is given by 

i^'^ix^y) = ycos"^(T(x,ro))/ro (4.4.5) 

and the inverse to V' on 7^ x [ri, vr] (the inverse to (j4.2.6|) ) is given by 

■ip~^(x,y) = cos^"'^(r(x, ri)) + ^(vr — cos~"'^(r(x, ri))). (4.4.6) 

vr — ri 

We will now prove some bounds on the derivatives of 'tp~^. 
Lemma 4.4.5. Z^i^'^ G [-^'O]- 

Proof. Let (x,y) G 7^ x [0,ro], so tp^^ is given by (j4.4.5p . Then 



Di^p-\x,y) 
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where r = r(j;,ro). Using Lemma 14.4.21 the numerator equates to — y/4, having range 
[—1/2,0]. Similarly the denominator has range [3/2, \/3], giving Diip^^ G [—1/3,0]. 
For {x,y) £ TZ x [ro,7r] the function ip^^ is given by (j4.4.6p and so 



(vr — ri)\/l — ' 



where r = r(x, ri). We find that the numerator is in the range [3(\/7-7r)/4,0], whereas 
the denominator is in [ll(7r - V7)/12,tt - V7]. Thus Diip"^ G [-9/11,0]. 
Lastly if {x,y) £ TZ x TZ then ip'^ is given by (j4.4.4p and so 



Diij-\x,y) 



X 2 



72' 



where r = T{x,y). The numerator equates to — 1/4 + (4 — y^)/4x^, which is increasing 
in X (because ^ 4) and decreasing in y. It has the range [-7/16,-1/4]. Using 
Lemma [4.4.2l and the fact that r is increasing in r, the denominator is in [3/4, 1], giving 
Di^/^-i G [-7/12,-1/4]. □ 
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Lemma 4.4.6. D21P ^ G 
Proof. Let £ TZ x [0, tq], then 

D2ip~^{x,y) =cos"V/ro, 

where r = r(2;,ro). From the proof of Lemma 14.4.41 we have 021/^^^ £ [7/12,3/4]. 
Next, let (x,y) £ TZ x [ri,7r], then 



D2^^J ^{x,y) 



IT — cos ^ T 



vr — ri 



with T = T{x,ri). Comparing with the proof of Lemma 14.4.41 we have D21P ^ £ [3,4]. 
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Finally let {x,y) G then 



Hx,y) = - — . 
2x\/l — 

where r = T{x,y). In the proof of Lemma 14.4.51 we established that \/l — G [3/4, 1] 
and so the denominator is in the range [3,2^/7]. Thus we have D2^p-^ G [V7/7,V7/3]. 

□ 

With the bounds established in Lemmas 14.4.31 through 14.4.61 it is a simple matter 
to check that the conditions (j4.4.3p will always be satisfied, and thus the proposition is 
proved. 

4.5 The Bernoulli property 

In this section we improve our estimate on the direction of unstable manifolds in the 
following sense. Let z = {x,y) £ R. We show that E^{z) C C C TzR and so 7"(z) is 
aligned within the cone C (we will say precisely what we mean by this in a moment, 
although the reader probably has an intuitive idea). Using this fact we are able to 
deduce that the strong form of the manifold intersection property, condition (j4.3.2p . is 
satisfied. 

4.5.1 Orientation of the unstable subspace 

Let uj G M^"^(E) be 'typical' in the sense of Section which is to say that ■j'^{uj) 
exists, its length is well defined and for any n G N the length of Q'^{j'^{uj)) is well 
defined also. Recall that in this case the latter grows exponentially with n. Let I be an 
interval and a : / — > MJ^^{A) the smooth curve whose trace is j^(oj). 

Let z = ^(uj) give such a point in the new coordinates, and consider j'^{z) = 

^(7"(^))- 
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Definition (Alignment of local invariant manifold). We say that "y'^{z) is aligned within 
the cone C , or that it has orientation in the cone C , if and only if the derivative 
D^Dao G C C T^R. 

Suppose that 7" (2:) is aligned within the cone C. It follows immediately from Propo- 
sition [53?T] that 0^(7"(u;)) is aligned within the cone C also for any n G N. This is 
enough for us to prove the Bernoulli property as below. 

Unfortunately, verifying that our assumption holds will require o ne last digression. 



For a full-measure set w G A, denoting uj = Ml^^{w) then the results of lWoitkowski (19801 ) 
tell us that 'J^{uj) has orientation in the cone U. It does not necessarily follow from our 
coordinate transformation that if z = ^(to) then 7^(2;) has orientation in the cone C. 
To overcome this problem we formulate the following proposition: 

Proposition 4.5.1. Let X be a compact metric space of dimension 2 and T : X ^ X 
a (non-uniformly) hyperbolic transformation, preserving a measure fj, on X . Let Y C X 
be a subset to which fx-a.e. trajectory returns infinitely many times (i.e. {T"(x)}„eNny 
is infinite) and for which the 'first-return map' Ty :Y is uniformly hyperbolic. 

Suppose that T is differentiable fi-a.e. and define cones K{x) C TxX of the form 

K{x) = {(r/,C) : ki{x) < C/r] < A;2(x)} , 

where ki,k2 G [—00, 00]. If DTx {K{x)) C K{T{x)) and the boundary of K{x) is mapped 
to the interior of K{T{x)), then £'"(x) must lie in the interior of K{x). 

Proof. In the tangent space T^X define unit vectors in the stable and unstable subspaces 

[si[x),S2{x)) G E\x) and {ui[x),U2{x)) e E'^{x) 

respectively. Fix x G X and let (r^O) Co) S K{x) C TxX. There are unique (and, without 
loss of generality, non-negative) real constants ao and [3q such that 
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(%,Co) = ao{si{x),S2{x)) + /3o{ui{x),U2{x)). 



Similarly, if we fix n G N and consider (r7„,Cn) = DT^{rio,Co) e K{T''{x)) C Ttu^^-^X 
then there are unique and non- negative an,(3n such that 

iVnXn) = Msi{T^{x)),S2{T^ix)) + Pn{uiiT^ix)),U2{T''{x))). 

The uniform hyperbolicity of the return map to Y C X allows us to estimate the 
magnitudes of an-,Pn- Given m G N and x G F there exists some n G N such that 
r™(x) = r"(x). Under such circumstances we have 

an^A-'^ao and /3„ ^ A"*/?o, 

where A > 1 is a constant independent of m, n or x. It must be the case that m — ^ oo 
as n — > oo (though we say nothing about the relative rates of divergence) and so as 
n ^ oo we see that and /?„ oo. 

K{x) cannot be contained in E^{x), so there is {r]o,(o) G K{x) for which /3o / 0. It 
follows that 

Cn _ U2{T^{x)) + (a„//3„)s2(T"(x)) 
r/„ ui{T^{x)) + (a„//3„)si(r'^(x)) ' 

where aJPn ^ X'^'^ao/Po ^ as n — oo. Consequently 

|Cn/r?n-«2(T"(x))/ni(r"(x))| 

tends to zero and, because {r]n,Cn) is in the interior of K{T'^{x)), for all sufficiently 
large n, so is £'"(r"(x)). The fact that x was arbitrary completes the proof. □ 
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4.5.2 The Bernoulli property 

It remains to show that the divergence and orientation of manifolds we have estabhshed 
give sufficient conditions for the Bernoulh property. This amounts to showing that the 
strong form of the condition (|4.3.2p is satisfied. 

Just as we have estabhshed that for 'almost every' z £ R the unstable manifold 
7"(z) grows exponentially when iterated with Hs and its orientation remains within 
the cone C, so it can be shown that for 'almost every' z' €z R the stable manifold 7^(2;) 
grows exponentially when iterated with H^^ and its orientation remains within the cone 
C. 

We need to develop a little terminology. In particular we will introduce a covering 
space for the manifold R, onto which we will lift our local invariant manifolds and deduce 
intersections. Any intersection of the lifted local manifolds must imply an intersection 
of the local manifolds themselves. We construct the covering space in two stages. 

Let —R C denote those points (x, y) C such that (— x, —y) G R. We notice 
that Rn —R = and define a manifold Ri = RU —R. Figure [491 illustrates Let 
p' : Ri ^ R he given by (x,y) 1-^ {x,y) if (x,y) G R and (x,y) 1— > (— x, — y) otherwise. 
Then is a covering space (a double cover, in fact) of R. The derivative of p' and 

its possible inverses preserve the cones C and C. 

Recall that {u,v) give Cartesian coordinates in the plane M^. We define 

R2 = {{u, v) : rQ ^ \u — 2n7r|, \v — 2rmi\ ^ ri, for some m,n G 7^} . 

A portion of R2 is illustrated in Figure H. 101 Let p" : R2 ^ Ri be the projection which 
takes each coordinate modulo S^. Then {R2,p") gives a covering space for Ri. The 
derivative of p" and its (locally defined) inverses preserve C and C. Thus {R2,p = p'op") 
is a covering space for R, and Dp, Dp^^ preserve C and C. 

Lemma 4.5.2 (Lifting lemma). Let z £ R and let p be a curve into R, that is, a 
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-IT —ri — ro 



ro ri TT 



Figure 4.9: The manifold Ri = RU R' CT'^. Together with the map p' : Ri 
gives a covering space for R. 



R this 



continuous map from some interval containing into R. Furthermore let p{0) = z. 
Then for any z" G R2 lying over z (i.e. p{z") = z) there is a unique curve p" into R2 
lying over p (i.e. p{p") = p)- 

Most topofog y texts contain a proof of this result; we recommend the book of 
Armstrong (19831 ) . We remark that what we have called curves are often called paths 
in the topology literature. 

We are ready to prove the main result. 



Katok et al. (19861 ). It suffices 



Proof of Theorem \1.3.2l Recall Theorem 12.2.51 due to 
to show that the strong form of the manifold intersection property (j4.3.2p is satisfied. 
Expressed in the new coordinates, the condition states that for p-a.e. pair w,w' £ A 
(writing z = o M^^{w) and z' = ^ o M^^{w')) and for all sufficiently large integers 
m and n then 

/7"(7"(^))n/7-'"(7'(z)) /0. (4.5.1) 

It has been established that if" (7^(2)) has orientation in C and that so do any of its 
lifts to i?2. By choosing n large enough a lift may be as long as we would like. Moreover 
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Figure 4.10: A portion of the manifold i?2 C M^. Together with the map p : R2 ^ R 
this gives a covering space for R. In red is a typical piece of the image of a local unstable 
manifold for some z £ R. The gradient at all times is in C. Analogously we show a 
typical piece of the pre-image of a local stable manifold for z' G R. This has gradient 
in C. If they are sufficiently long then they must intersect. 
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the lifting lemma tells us that any of these lifts is a line with direction 'bottom-left' to 
'top-right', as we have illustrated the manifold in Figure H.lOl 

Whenever j^{z) crosses an intersection region (marked -|- or — ) horizontally (re- 
spectively, vertically) then the next iteration with F (respec. G) stretches its image 
vertically (horizontally) adding at least 2ir to its height (width). This shows that the 
pathological case of H^{'^'^{z)) being oriented in C but essentially horizontal or vertical, 
cannot occur. 

We are left with H^'{"f^{z)) having arbitrary length and stretching diagonally as 
previously described. Such a piece is illustrated in red in Figure 14.101 Analogous 
arguments show that H~'^{"f^{z')) must be as illustrated in blue and an intersection is 
inevitable. □ 
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5 The Bernoulli property for a 
linked-twist map on the 
two-sphere 



In this chapter we prove Theorem 11.3.31 which says that the hnked- twist map on i? C 
defined in Section [1.3.31 has the Bernouhi property. 

In Section [5. II we introduce a generalised linked-twist map H : R ^ R on a manifold 
i? C T^. We append the word 'generalised' because the new map does not fit the 
definition of an abstract linked-twist map given in Section 11.2.21 consisting as it does of 
twist maps defined on four annuli rather than two. Each of the twists is linear. 

We are able t o prove, by application of the techniques intro duced by 



and amended by 



Woitkowski (IdSOD 



Katok et al. (IdSGi ). that 



Sturman et al. (20061 ) in light of the work of [ 
the new map has the Bernoulli property. In Section 15.21 we deal with some technical 
issues concerning the nature of points at which H is non-differentiable and conclude 
that Lyapunov exponents exist Lebesgue-almost everywhere. We then describe in some 
detail the return of points to a certain region S C R before using this description to 
show that these Lyapunov exponents are non-zero. 

In Section 15.31 we deal with the 'global' aspects of the argument. We consider 
the orientation of local invariant manifolds and give a rigorous original proof of their 
growth. Such a result is lacking in the literature and is similar to the analogous result 
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from Chapter m with it we are able to prove the Bernoulh property for H. 

In Section 15.41 we show that H is semi-conjugate to the linked- twist map on §^ 
and from this c onclude that the latter is also Bernoulli. Our proof relies on a result of 
Ornstein (197ll ). 

5.1 A generalised linked-twist map on the two-torus 

In this section we define a more general linked-twist map on the two-torus T^. It differs 
from the map introduced in Section 11.3.11 in that it is constructed by embedding not 
two but four cylinders into T'^. 

We will define the generalised linked-twist map directly on the two-torus; we take it 
to be elementary that the manifold we construct may be obtained by embedding four 
cylinders into and do not give a formal proof. 

5.1.1 Definition of the map 

Recall our construction in Section 11.3.31 of the linked-twist map Q : A ^ A, where we 
described the circle as having a coordinate that is periodic with period AK and where 



Presently it will be most convenient to let this coordinate take values between —K and 
2>K, where of course these particular values are identified. Let < XQ^y^ < K he as 
previously defined. We define two 'horizontal' and two 'vertical' annuli on the two-torus: 



Denote by P = PqUPi the union of the horizontal annuli and hy Q = Qo^Qi the union 
of the vertical annuli. We let R = PUQ and Shi = Ph^Qi for each pair h,i £ {0, 1}. Let 




(5.1.1) 



Po = {{x,y) : X G S\ \y\ ^ yo}, 
Qo = {{x,y) : \x\ ^ xo,y G 



Pi = {(x,y) : X G SMy - 2KK yo}, 
Qi = {{x,y) :\x-2K\ ^xo,yeS^}. 
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S denote the union of the four regions Shi. Figure [STT] illustrates the manifold i? C T^. 

Qo Qi 



2K 




















Su 
















Sm 




Soi 















Pi 



-Xq Xq 2K 
Figure 5.1: The manifold i? C is shaded. 

We define a twist map on each of the four annuli and group these into two 'families'; 
those maps defined on horizontal annuli are in the first family and those defined on 
vertical annuli are in the second. Both maps in a given family act simultaneously and 
the linked-twist map will be the composition of the two families. 

Let / : ^ (see Figure 15. 2p and 5 : §^ — > be given by 



f{y) = { 



iK{y + yo)/{2yo) if y G [-yo, yo], 
4K{y + yo- 2K)/{2yo) if y G [2K - yo, 2K + yo], 
otherwise, 



9[x) 



4:K{X + Xq)/ (2xo) if 2; e [-Xq, Xq], 

AK{x + X0- 2K)/{2xo) if x G [2K - xo,2K + xq] 



otherwise. 

A family of horizontal twist maps F : — > is given by 



F{x,y) = {x + f{y),y) 
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4K 



2K 








-K -2/0 Vo 2K-yQ 2K+yQ 



Figure 5.2: The generalised twist function / : S is continuous everywhere and 

difFerentiable except at the four points ityo and 2K it i/q. The function (/ : ^ is 
similar. 

and a family of vertical twist maps G : — > is given by 



G{x,y) = {x,y + g{x)). 



We say that the map H : R ^ R given by the composition H = G o F is a gener- 
alised linked-twist map on R. By analogy with the toral linked-twist map introduced in 
Section 11.3.11 it is easy to see that H preserves the Lebesgue measure on R, which we 
will denote by fi. In Sections 15.21 and 15.31 we prove the following: 

Theorem 5.1.1. The map H : R ^ R is metrically isomorphic to a Bernoulli shift. 



5.2 Non-zero Lyapunov exponents for H 

A major step toward the proof that H has the Bernoulli property is to show that its 
Lyapunov exponents are non-zero. The present section is divided into three parts. In 
Section [5.2. II we g ive some results of a technical nature which satisfy certain hypotheses 
in the theorem of lKatok et al. (19861 ) (Theorem 12. 2. 5p . This allows us to conclude that 
Lyapunov exponents exist ^-almost everywhere. In Section [5]2]2] we describe the return 
of a 'typical' trajectory to the region S; this behaviour will be crucial in determining 
hyperbolicity. Finally in Section [5.2.31 we describe the behaviour of the differential DH 
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acting on certain tangent cones, from which the result fohows. 



5.2.1 Technical details 

Our first consideration is to determine the nature of those points in R at which H and 
its iterates are non-differentiable. For this we will need to establish some notation for 
boundaries of the annuli. The set 

dPo = {{x,y) :xeS^,ye {yo,yi}} 

denotes the boundary of annulus Pq, with the boundaries of Pi, Qq and Qi being 
defined similarly. The boundaries of P and Q are denoted by dP = dPo U dPi and 
dQ = dQo U dQi respectively. 

The twist function / is differentiable provided that y,y — 2K ^ ityo and similarly 
g is differentiable provided that x,x — 2K ^ ±xo. It follows that F is differentiable 
on T^\9P and that G is differentiable on T'^\dQ. Consequently H is differentiable on 
T^\so, where 

so = dPUF~^{dQ). 

It may be useful for the reader to briefly review Secti on |2.2.3l at the present moment, 



where we described the theorem of 



Katok et al. (19861 ) and the notation necessitated 



by its statement. We remark that the function / from that section corresponds to 
the function H from this section; the function / in this section has no analogue in 
Section [2X3l 

We define a number of full-measure subsets of R and hope that Table 5.1 is of use 
to the reader in keeping track of these. 

Let p he a metric on R, so that {R, p) is a complete metric space. Notice that 



^This notation is necessarily different from that in Chapter |3] where both our map and our ambitions 
differed. 
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Set 



R 
V 
N 
J 



Description 



Union of the four annuli 
Open subset; a Riemannian manifold 
Open subset on which H is defined 
Intersection of ah images and pre-images of N 



Table 5.1: Full measure subsets of the complete metric space R. As listed, each set 
contains those below it. 

V = R\sq is a smooth manifold. V has full /n-measure and is open. Let N = F\iJ~^(so); 
N also has full //-measure and is open. 

Proposition 5.2.1. H\jq : N ^ V is a smooth map with singularities. 

Proof. N is open, so given any z E N the exponential map exp^ is an injective map 
from some neighbourhood of G T^T^ to some neighbourhood U{z) C N of z. The 
only restricting factors on the size that this neighbourhood may take are the finite size 
of N itself and the distance from a given point to the set sing(iJ) = V\N. Furthermore 
if z e N then z ^ dP and F{z) ^ dQ so H is smooth and injective. □ 

In our next lemma it will be convenient to consider the Lebesgue measure on D i? 
and we denote this by iJ.f2. For //t2 -measurable sets C T^, /x (the Lebesgue measure 
on R) is the conditional measure such that 



l^iK) = fij2iK\R) = 



Hf2 {K n R) 

HY2{R) 



It follows that for measurable sets K' (Z R(ZT'^ we have 



(5.2.1) 



We use the notation Bj2{z,e) to denote the open e-neighb our hood of z G and 
we use B{z,s) to denote the intersection Bj2{z,s) fl R. For a non-empty set K C T'^ 
we let 
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Bj2{K,e) = U Bj2{z,e) 

and similarly define B{K,e). 

Lemma 5.2.2. H\fyf satisfies the condition i2.2.10\) . 

We also refer to this as the condition (KSl). 

Proof. We have smg{H) = H-^{so) = H'^ {dP U F-^{dQ)) . In T^, an e-neighbour- 
hood of dP consists of four strips (rectangles) each of width 2e and length 1, thus 
fj.j2{Bj2{dP,e)) = 8e. B{dP,e) is a subset of Bj2{dP,e) and so its //Tr2 -measure is at 
most 8e. We are interested in its /i-measure and, by (j5.2.ip above, these are related by 

fij2{B{dP,e)) = fi{B{dP,e))fij2{R), 

all of which gives ^{B{dP,e)) ^ 8e/fij2{R). 

F~^{dQ) consists of four piecewise-straight lines, each of length bounded above by 
some integer A^. Thus a similar argument applies here. Taking the pre-image with 
respect to H will 'stretch' these strips somewhat, but similar bounds will still apply so 
that (KSl) will be satisfied with some constant c G M and with a = 1. □ 

Denote by DF, DG and DH the derivatives of F, G and H respectively. We use 
subscripts to denote the point in R at which the derivative is evaluated, so for example 
at each z ^ N we have the familiar chain rule 

DH, = DGf(^,)DF,. 

Let J = Pli^-oo H^{N); we prove that this has full ^-measure: 
Lemma 5.2.3. fi{J) = 1. 
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Proof. We use the ff-invariance of /i and the countable additivity property. From 
elementary set theory we have J = R\{R\J} and so 



MJ) =KR)-KR\J) 

^KR)-Er=-oof^{RW{N)) 
= KR)- 

The fact that J C R completes the proof. □ 
Lemma 5.2.4. Lyapunov exponents for H\]\f exist for fi-a.e. z € R and every non-zero 

V e r^T2. 

Proof. Let z = (x, y) G J. The derivatives of F and G are given by 



1 f'{y) \ / 1 

DF=\ and DG = 

1/ \ g'{x) 1 



respectively. Each takes only finitely many values on R meaning that there are upper 
and lower bounds on its eigenvalues. Thus 

/ \og^ \\DH^^\\dn= [ log^ \\DH^^\\dn < oo 
Jr Jj 

(the norms are operator norms), showing that H satisfies (j2.2.12|) . The result follows 
immediately from the multiplicative ergodic theorem, Theorem 12.2.41 □ 

Lemma 5.2.5. H\n satisfies the condition \2.2.11\) . 

We also refer to this as the condition (KS2). 

Proof. Because our twist functions / and g are piecewise-linear, their derivatives (where 
they exist) are constant functions. Consequently all second derivatives will be zero and 
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condition (KS2) is satisfied with b = I and any C2 > 0. □ 

Lemmas 15.2.21 15.2.41 and 15.2.51 and part (a) of Theorem 12.2.51 yield two conclusions. 
First, Lyapunov exponents exist for ^u-a.e. z € i? and for every non-zero tangent vector 
at z. Second, for those z at which we have a positive exponent, the local unstable 
manifold 'y^{z) exists and is of the form (|2.2.9p . An analogous statement holds for 
negative exponents and the local stable manifold j^{z). 

Throughout this section we ha ve referred to the restriction H\fyf in keeping with 
the notation of 



Katok et al. (19861 ). We have seen that N has full /x-measure and so 
statements that are true for a.e. z £ N are equally true for a.e. z £ R. From here on we 
will speak only of H rather than its restriction to A^, this leading to cleaner expressions 
whilst losing none of the accuracy. 

5.2.2 Return of trajectories to the region S 

We will need to describe in some detail the trajectory of a 'typical' point z £ R (we will 
say precisely which points are 'typical' soon; for now let us just say that the set of such 
points will have full measure). In particular we need to analyse first whether, and if so 
how often, such a point returns to the union of intersection regions S. Our main tool 
for doing this will be the return map to the set S, which is defined as follows: 

Definition (First (and n^^) return map to S for z). For z £ R, let n he the smallest 
(strictly positive) integer such that H^(z) £ S, if such an n exists. If it does then the 
map 

Hs{z) = H^{z) 

will he called the first return map to S for z, often abbreviated to just first return map 
where it is clear for which z the map is defined. Similarly 

Hl{z) = Hs{H^{z))oHs{z) 
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will be called the second return map and the general case follows inductively. 



In the remainder of this section we describe the notation and formahsm introduced 



by IWoitkowski (19801 ) in his study of toral hnked- twist maps defined on two annuli. It 
win enable us to give the description we require. Let p,q : R ^ [0, oo] be defined as 
fohows: 



p{z) 



if z e Q\P, 

pen ifzGP, ^ Q for < n < p and G Q, 

oo if z e P and F"(z) ^ Q for all n G N, 



if z G P\Q, 
g G N iiz£Q, ^"(z) ^ P for < n < g and Ciz) G P, 
oo if z £ Q and G"(z) ^ P for all n G N. 

We give a brief discussion. Let z G P, so p{z) ^ and assume also that p{z) ^ oo. 
Clearly F^{z) G P for all natural numbers n, and p{z) is the first such natural number 
for which F^{z) G Q, i.e. the first time the forward trajectory under F lands in S. This 
is of interest because for < n < p{z) we have F'^{z) ^ Q and so G{F^{z)) = F^{z). 

To phrase this another way, starting with z as above, if we follow its orbit under 
H = G o F then the p{zY^ iterate is the first for which the G component of H is not 
the identity. A similar property holds for the function q. 



Wojtkowski (1980l ) outlined an argument that the trajectory of fi-a.e. z £ R returns 
to S infinitely many times. The proof is quite elementary and we have discussed an 
analogous result in Section 14.11 We denote by P^ C R the full measure (and clearly 
invariant) set of points for which this occurs. 

Define two sequences of functions Pi '■ Rg ^ oo), pi : Rg ^ [1, oo) for i = 2, 3, 4, ... 
and Qi : Rg — > [1, oo) for i = 1, 2, 3, ... as follows^ 



^We will often neglect to express the dependence of Pi,qi on the point z for the benefit of cleaner 
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Pi+i{z) =p{GiiFPi ■■■Gi^FP^{z)), qi+i{z) = qiFH+^GUFPi ■ ■ ■ G1^FP^{z)). 



Given z G R§ the above scheme defines inductively pi and qi for all natural numbers i. 
It is clear that for z e Rg one has e S; for z e R§ D S then ^^1(2;) G S 

also. 

Finally define functions mi : Rg — ^ [0, 00) and mi : Rg — > [1, 00) for i = 2, 3, 4, ... by 

i 

n=l 

The description of the first-return map that we required is given by the following easy 
lemma: 

Lemma 5.2.6. For z G Rg, the first return map to S for z is given by Hs = H'^'^^^\ 
In general, for n G N, the n*'* return map to S is given by Hg = ffi^niz) _ 

Proof. We have 

H""! = (G o F) o • • • o (G o F) o (G o F) o • • • o (G o F) 

= {G o id) o ■ ■ ■ o {G o id) o {G o F) o {ido F) o ■ ■ ■ o {ido F) 
qi-1 Pi 

i.e. H"^^ = G^i o F^'i . The first result follows easily from the definitions of pi and q^ and 

the general case by induction on n. □ 

5.2.3 Lyapunov exponents are non-zero 

In this section we use the first return map to show that Lyapunov exponents are non- 
zero almost everywhere. Our first step is to describe the derivative of the first return 
map. 



notation, where this will not lead to confusion; for the same reason we typically neglect the composition 
symbol 'o'. 
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Let z G n J and denote z' = FP'^^'^\z), then the derivative of the first-return 
map is given by DHs = DH^^ = DG'i,] DFi^ , where pi,qi and rrii are to be evaluated 
at z. Similarly for n G N the derivative of the n*'^ return map is given by 

where Zq = z, z'q = z' and in general = H"^^{z), z\ = for i € N. 

Let (u, t;) = (dx,dy) give coordinates in the tangent space T^T^, which we identify 
with M^. We define the cone 

C = {{u,v) G : ui; > 0} C T^T^ 

that is the open first and third quadrants of the plane, and introduce the standard 
Euclidean norm || • || : ^ [0, oo) on T^T^. 

We show that for /x-a.e. z ^ R and for each w G C C TjT^ the Lyapunov exponent 

X'^{z,w) = lim — log 

n— »oo fi 

is positive and thus there is a local unstable manifold 7^(2;) of the form (12.2.9p . We do 
not do so explicitly, but the results are easily re-formulated to show that if 

C = {{u,v) G : TO < 0} C T^T^ 

then for fi-a.e. z € R and for each w £ C the Lyapunov exponent 

X^iz,w)= lim —log\\DH^w\\ 

n—*~oo n 

is less than zero, and so a local stable manifold 7'^(-z) exists also. 
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The first step is to establish some results concerning the growth of a vector w € C 
when acted upon by the derivative of H or of Hs- Similar results were proven in 



Woitkowski (198d l. The result illustrates a key feature of the dynamics: the return 
map is uniformly hyperbolic. 

We use the convention that N denotes the positive integers; in particular ^ N. 

Lemma 5.2.7. Let z € D J and let w e C C T^T^. Let n, r G N with r ^ n. Then 

(a) DH^w eC C rf/n(^)T2, 

(b) \\DHlw\\ ^ \\DHl^w\\. 

Lf additionally z £ S then for a constant A > 1, which is independent of z, we have 



(c) 



> A"||w||. 



We remark that part (a) of the lemma says that the cone C is invariant under 
the derivative of and that part (b) is a 'non-shrinking' condition, saying that the 
norm of vectors in C does not decrease when mapped by DH. Part (c) expresses the 
uniform growth resulting from the derivative of the return map. The extra condition 
that z £ S will not limit our ability to make use of this lemma, because we are interested 
in cumulative growth along a trajectory and we have already established that /i-a.e. 
trajectory returns infinitely many times to S. 

Proof Let z = {x,y) £ Ra H J and let z' = (x', y') = F{z). We have 




where we have introduced the notation a = f'{y) and /? = g'{x'). At most one of a, /3 
may be zero and the other(s) positive, because a = requires (x, y) ^ P and (3 = 
requires F{x, y) ^ Q, which cannot occur in succession. 
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Now let w = {u,v) £ C C T^T^; straight-forward calculation shows that DHzW S C 
and that the Euclidean norm of this vector is at least that of w. Induction on n yields 
the result (a) and that ^ \\w\\ for n G N. The result (b) follows easily. 

Now assume that z €z S and redefine z' = {x',y') = FP^^^\z). (The functions 
Pi, Qi, rrii are always to be evaluated at z.) Notice that (the redefined) a, f3 are now both 
strictly positive. We have 



1 pia 
qi(5 l+piaqif5 



where a = f'{y) > is a positive constant, as is /3 = g'{x') > 0. Let k = min{pia, qi(3} > 
and A = VTTk^ > 1. Evaluating \\DH^^ wW gives 



u2(l + qf/3^) + v^ipja^ + (1 + PiaqiPy) + uv{2pia + 2gi/3(l + piaqi/3)) 



for which a lower bound is (u^ + v'^){l + k^) = A||t(;||. This proves part (c) in the case 
n = 1. The general case follows by induction on n. □ 

The previous lemma shows that the tangent cone C is uniformly expanded by DHs, 
but not necessarily by DH. In order to demonstrate the non-vanishing of Lyapunov 
exponents we must consider the frequency with which a typical trajectory returns to 
the intersection region. To that end we use the following result of Burton and Easton: 



Lemma 5.2.8 (iBurton and Easton (19801 )). Let X be a compact metric space, fi a 



Borel measure on X and T : X ^ X a fi-preserving homeomorphism. Suppose Y C X 
is measurable. For /x-a. e. x £ X such that T^{x) G Y for some n E N, that forward 
orbit will return to Y with positive frequency in n. 

The result we want is an easy corollary of Lemma 15.2.81 and the fact that /u-a.e. 
z £ R enters S: 
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Corollary 5.2.9. For fi-a.e. z £ R, the forward orbit 0^{z) = {H"-{z) : n G N} lands 
in S with positive frequency in n, i.e. there is a set Rs C R^ C R, li{Rs) = !> o-iT-d for 
each z G Rs the limit 



exists and is strictly positive. 

We remark that Lyapunov exponents (where they exist), being infinite-time hmits, 
are invariant along a given trajectory. That is, if x{^i exists for some z £ R and 
non-zero w G T^T^, and if n G N, then x{H^{z), DH^w) exists and is equal to x{z,w). 
The main result of this section is the following. 

Proposition 5.2.10. For fi-a.e. z £ R and for every tangent vector w £ C C T^T^ at 
z, the Lyapunov exponent x^iz,w) for the map H is positive. 

Proof. Let z £ RsDJ and let u; G C C T^T^. We denote z' = (z) and w' = DH^^w. 
Clearly z' £ RsCiJnS and by Lemma r5.2.7l we have w' £ C C T^/T^. By the invariance of 
Lyapunov exponents along a given trajectory it is equivalent to show that x^{z' , w') > 0. 

We claim that for any z' £ Rs f) J D S there is a positive constant and if n > 
then n > m^^^^, where r(z) > is a constant to be determined. This condition means 
that for the trajectory of z' , after some 'transition period' given by N, there is a positive 
lower bound on how frequently it hits S. 

Suppose for a moment that our claim is verified. Lemma 15.2.71 implies that, corre- 
sponding to each such return, the tangent cone C is expanded by a factor A > 1. The 
situation is as follows: 



6{z) = lim 



n 



w') = lim. 



■n- 



oo 



l\og\\DH^,w' 



^ lim, 



^ lim. 



■n- 



■n- 



oo 



oo 



ilog DH^-^w' 
ilogAL™J||u;'|| 



r log A > 
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where the second hne follows from our claim (we may, of course, assume n > in the 
limit n ^ oo) and from Lemma 15.2.71 the third from the same lemma and the fourth 
from elementary properties of the logarithm. 

It remains only to prove the claim. Fix z' £ Rs Ci J D S. By Corollary 15.2.91 there 
exists d{z') = \imn~>oon/mn{z') > 0. In particular, there exists > such that if 
n > N then n/mn{z') > 6{z')/2, or equivalently that 2n/6{z') > mn{z). If we make 
the substitution n' = 2n/5{z') then it follows that there exists N' = 2N/5{z') > such 
that for all n' > N' we have n' > m^n'S{z')/2} > which is the statement of our claim, with 
r{z') = 6{z')/2>0. □ 

We have remarked that it is completely analogous to show that for /x-a.e. z £ R 
and for each non-zero w £ C C T^T^ the Lyapunov exponents x~ (-^i u)) are strictly 
nega tive. R has dimension tw o, so by the elementary theory of Lyapunov exponents 
(see 



Barreira and Pesin (20021 )) there are at most two distinct Lyapunov exponents 
associated to z and we have shown them both to be different from zero. 

Part (b) of Theorem 12.2.51 implies that our system has an ergodic partition and that 
moreover, associated to /i-a.e. z G R, there is a local unstable manifold 7"(2;) and a 
local stable manifold "y'^{z). 



5.3 Global arguments 

We conclude the proof of the Bernoulli property by giving the 'global' aspects of the 
argument. We have established that H : R ^ R has an ergodic partition and it 
remains to show that there is just one positive measure component (i.e. a full-measure 
component) to this. We begin by describing the orientation of local invariant manifolds 
and deducing that their lengths diverge on iteration with H (or H~^, as is appropriate), 
then show that this behaviour is enough for us to conclude that H is Bernoulli. 
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5.3.1 Orientation of local invariant manifolds 

We discuss the length of local unstable manifolds and of their iterates under H. The 
situation is analogous to Section 14.31 in which we looked at local unstable manifolds 
for the planar linked-twist map. Here our task will be simplified by properties of the 
tangent cones C and C. 

For ^-a.e. z £ R the local unstable manifold J^{z) exists and is of the form (j2.2.9p . 
Thus there is an open interval / C M and a diffeomorphism a : / ^ 7^(2) C R. At z, the 
local unstable manifold is tangent to E^{z) C T^. It is clear (although Proposition 14.5.1] 
may be used for a rigorous proof) that E'^{z) C C. 

In order to define the length of iterates of the local unstable manifold we need the 
following lemma: 

Lemma 5.3.1. For fi-a.e. z £ R and for any n £ N, the set of z' £ 'y^{z) at which 
DH^ does not exist is at most countable. 

Proof. The proof is a consequence of the orientation of 7^(2;) and of the constituent line 
segments of the set of non-differentiable points for DH^. The former, we have argued, 
lies in C. We claim that the latter is given by the set 



We justified this previously in the case n = 1 and the general case follows by induction 
on n. 

Figure [531 illustrates the set sq. We did not construct local stable manifolds explic- 
itly, but had we done so we would have formulated a result completely analogous to 
Lemma l5.2.7l showing in particular that the cone C = {{u,v) : < 0} is preserved by 
DH^^. By referring to the proof of Lemma 15.2.71 it is straight-forward to see that the 
statement still holds if we replace C with its closure. 



n 




(5.3.1) 



i=0 
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-Xq Xq 



Figure 5.3: The manifold i? C showing those points (in red) at which H is not 
differentiable. We illustrate only those points in the interior of the annuli. The map is 
not differentiable on the boundary either. 

We conclude as follows. The non-differentiable set for H consists of finitely many 
line- segments, each having orientation within the closure of the cone C. These orien- 
tations are preserved by DH^^ so the statement holds true for the non-differentiable 
set for H"^, with n some positive integer. Because 7"(-z) has transversal orientation, all 
intersections 



are transversal. It follows from Lemma 14.3.31 that there are at most countably many 
such intersections. □ 



Our proof of the growth of local unstable manifolds is now very similar to that of 
Theorem 14.3.11 



Theorem 5.3.2. Let n G N and let X > 1 be the constant given by Lemma 5.2.7 For 
jjL-a.e. z G R the length of Hg{z) is defined and 



length (H'^ij^'iz))) ^ X"" length ij"" (z)) . 



(5.3.2) 
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Proof. For //-a.e. z £ R there is an interval I C M and a diffeomorphism a : I ^ 
7*'(z) C R. Lemma l5.3.1l says that for a.e. such z and for any n G N the interval I has 
a countable (at most) decomposition 

TV 

I = {io,k) U [J [ih,ih+i) 

h=l 

for some S N U {00} , and the restriction of Hg o a to any interval {i^, ih+i) in this 
decomposition is a smooth curve. The length is given by our definition above and we 
have 

length (i?g(7"(z))) = length(Fg o a) 

= \-j,\\Dat\\dt 

= A"length(a) 

= A"length (7"(z)) . 

□ 

We have deduced exponential growth of local invariant manifolds with respect to 
the first return map Hs, and that the orientation of these manifolds is restricted to the 
cone C. The Bernoulli property follows exactly as in Section [4. 5. 2[ This completes the 
proof of Theorem I5.1.1[ 

5.4 Proof of the main result 

In this final section we complete the proof of Theorem 11.3.31 as follows. The reader 
should recall the definition of an embedding given in Section 11.2.11 We show in Sec- 
tion 15.4.11 that E is an embedding of the cylinders C and C into §^ . In Section 15.4.21 
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we demonstrate that Q and H are semi-conjugate; this was after all the purpose of 
introd ucing the map H . We complete the proof in Section 15.4.31 by describing a result 



due to lOrnstein (19711 ) and proving some accompanying results of a technical nature. 

Before we begin we must caution the reader as to our notation. In previous sections 
(u, v) have been reserved for points in a tangent space. From here on (u, v, w) shall 
represent Cartesian coordinates in M^. Similarly whereas C has denoted until now a 
particular tangent cone, from now on it will denote the cylinder in which is (it will 
be proven) embedded into S^. 

5.4.1 E is an embedding 

We show that : C — > P C §^ is an embedding as defined above. By similar arguments 
one can show that E o N : C ^ Q C S'^ is also an embedding. We begin by reviewing 
some properties of the Jacobi elliptic functions which we will require. 



Proposition 5.4.1 (iMeyer (20011 )). Let sn, cn, dn : R ^ M denote the Jacobi elliptic 
functions, where the parameter k is fixed at \/2/2. Let K he as defined by i5.Ll\) . 

The functions sn and cn are periodic with period AK and dn is periodic with period 
2K. For any t gM. we have the relationships 

sn{2K -t) =sn{t), cn{2K - t) = -cn{t) and dn{2K - t) = dn{t); 

we have the addition formulae 



cn^(t) + sn^(t) = 1 and dn^{t) + ^sn^(t) = 1; 



the derivatives are given by 



— sn(t) = cn(t)dn(t), — cn(t) = — dn(t)sn(t) and —dn{t) = — -sn(t)cn(t); 
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sn is a homeomorphism of (—K,K), whereas cn is a homeomorphisni of {0,2K); and 
finally sn is positive on {0,2K), negative on {2K,AK), cn is positive on {—K,K), neg- 
ative on {K, 3K) and dn is positive for all t G M. 

Recall that it is convenient to consider C as a subset of . Our first result concerns 
the image of C with respect to E. 

Lemma 5.4.2. EiT"^) C S^. 

Proof. Let G and let {u,v,w) = E{x,y). Using Proposition 15 .4. l] we find that 

+v'^ +uP' = sv?{x)dT?{y) + cn^(x)cn^(y) + dn^(x)sn^(i/) 

= sn2(x) (1 - \sii^{y)) + (l - sn2(x)) (l - sii^{y)) 

+ (1 - |sn2(x)) sv?{y) 
= 1. 

□ 

Lemma 5.4.3. C C and c M? are smooth manifolds of dimension 2. 

Proof. To be precise we should say that C is a smooth manifold with boundary., although 
this introduces extra technical considerations and requires further definitions and we 
are not in fact interested in the behaviour of the twist map on the boundary, which 
consists only of fixed points. Thus we show that the interior of C is a smooth manifold 
of dimension 2. We justify this by observing that if we denote the boundary by dC C 
then it is clear that ^{C\dC) = n{C). 
Let < e < K and define subsets of 

Ui = {{x,y) : X G {-e,3K),y G {-yo,yo)} 

and 
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U2 = {{x, y):xe {K, AK + e),ye {-yo, yo)} . 



Let (pi : Ui —>■ C and (j)2 U2 —>■ C he given by 

<t>i{x,y) = (t)2{x,y) = {x mod S\y). 



It is not difficult to cfieck that Ui=i 4>i{Ui) covers tlie interior of C and that the other 
conditions given in the definition of a smooth manifold of dimension 2 are satisfied. 
We now turn to S^. Let 1/ C be the open unit ball in M^, i.e. 

y = {(x,y)GM2:x2 + /<!}. 

For t G {tt, V, w} define ■i/'± : F ^ by 

^± : {x,y) = -x^ - y'^,x,y^ , 



V'± : ix,y) = (^x,±Vl -x2 -y^,y^ , 



^± : {x,y) = (^x,i/,±vT^^x2^^) • 

Each pair {V'+(V^), V'-(^)} covers all but the great circle i = so the union of all 
six coordinate neighbourhoods certainly covers S^. The other properties are easy to 
establish. □ 

Lemma 5.4.4. E is a differentiable map. 

Proof. Let {x,y) G C, choose i £ {1,2} so that {x,y) G (piiUi) and choose t G {u,v,w} 
and either + or — so that E{x,y) G tp±{V)- Define 
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Ui = cl)r^oE-'oi;'^{V), 



(5.4.1) 



where E ^ denotes the pre- image (as we have yet to discuss the injectivity or otherwise 
of E). We must show that E{(j)i{U^)) C ^p±(y) and that 

h = {ip^^y^ o E o : Ul ^ (5.4.2) 

is differentiable at (f)~^{x,y). The former is merely a rearrangement of (j5.4.ip . From 
the definition of (pi and E it is immediate that they are differentiable (in the case 
of E we need Proposition I5.4.ip . The differentiability of h thus depends upon the 
differentiability of (V'±) ^ ■ Without loss of generality take t = u and choose +. We 
have (V'") ^ iu,v,w) = {v,w), which again is clearly differentiable. □ 

Lemma 5.4.5. E is an immersion. 

Proof. Let z = {x, y) G C. We must verify that DEz : T^C — > Tg^^)!^^ is injective. 
The Jacobian matrix for E will not be square and so we will not be able to check the 
injectivity of the differential DEz by simply computing its determinant (and showing 
that this is non-zero). 

However, the condition is equivalent to asking whether Dh^-i^^ ^-^ is invertible, where 
h is as defined in (j5.4.2p . for each admissible choice of i G {1, 2}, t G {u, v, w} and either 
+ or — . The Jacobian for h will be square so it will be simple to check its injectivity. 

Without loss of generality we take i = 1, t = u and take +. Then 

/i(0-i)(x,y) ={r±Y^ °E{x,y) 

= {tp±) ^ (sn(x)dn(y), cn(x)cn(y), dn(2;)sn(y)) 
= (cn(x)cn(y), dn(x)sn(2/)) . 

The Jacobian is given by 
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(— sn(x)dn(x)cn(y) — sn(y)dn(y)cn(x) 
— |sn(x)cn(x)sn(2/) cn{y)dn{y)dn{x) 

which has determinant 

— sn(x)dn(y) ^dn^(x)cn^(y) + -sn^(y)cn^(x)^ . (5.4.3) 

We require an expHcit expression for the domain of h: 

(PY^ o E-^ o tPKV) = cpl^oE-^ {{{u,v,w) G §2 . ^ > 0}) 
= (/>r'({(x,y) GC7:sn(x)dn(y) >0}) 
= 4>1^ {{{x,y) G C : X mod G (0,2/s:)}) 
= (0,2K) X [-2/0, yo], 

where in the third hne we have used Proposition 15.4.11 Now observe that sn(x) and 
dn(y) are strictly positive on the domain considered. The bracketed term in (|5.4.3p will 
be positive unless both dv? {x)ct\^ {y) and sv?{y)cn^{x) are zero. In particular we note 
that dT?{x)cn^{y) > for any x G M and y G {-K,K) D [-yo,yo]- □ 

Lemma 5.4.6. E is a homeomorphism. 

Proof. As in the proof of Lemma 15.4.51 it is simpler to work with the composition h 
defined by (j5.4.2p . From elementary results about compositions of functions we have 
that is a homeomorphism if and only if /i is a homeomorphism, for each admissible 
choice of i,t and it. Without loss of generality we again consider the particular case 
i = 1, t = u and +. We have 

h = ir+y' oEo^,:^-^ {E-\r+{v))) - v. 

Recall from the proof of Lemma 15.4.51 that in this case the domain of h is given by 
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(0,2i^) X [-yo,yo]- Let 



s = cn{x)cn{y) and t = dn(x)sn(?/) 

then using Proposition 15.4.1] we obtain 

= 2 ^1 - ^sn2(x)^ sn^iy) = (l + cn^ix)) sn^iy) (5.4.4) 

which, combined with = cn2(x) (l — sn2(?/)) gives 

2t2 + s2^cn2(x)+sn2(y). (5.4.5) 

We will solve these for 5 = sn'^{y) and C = cn^lx). Rearranging (I5.4.4P in terms of C 
and substituting into (|5.4.5p we obtain 

- {2f + - 1) C - = 

which we solve using the quadratic formula to get 



C = i (it^ + s2 _ 1 ± ^(2t2 + s2_l)2^452j . 



(5.4.6) 



It is clear that y (2*2 + ^2 - 1)^ + 4^2 ^ It"^ + s"^ - 1 with equahty if and only if s = 0. 
Clearly C = cn^lx) is non-negative and so in (j5.4.6p we must always take it to be +. 
We arrive at 



cn(x) = (^2*2 + s2 _ 1 + ^ (2*2 + s2 - 1)2 + 4s2^ (5.4.7) 

and again must determine the sign to be taken in (j5.4.7p . Notice (see Figure II. 7p that 
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cn(y) is positive for each y £ I = [—yo,yo]- Thus sgn(cn(j;)) = sgn(cn(x)cn(y)) = sgn(s) 
on the whole domain and cn{x) is well-defined by 

cn(x) = sgn(s)^^ ^2*2 + ^2 _ ^ _^ ^ (2f2 + s2 _ i)2 _^ 4^2^ _ (5 4 8) 

We need only observe (Proposition I5.4.T]) that cn is bijective on (0, 2K) to see that x is 
uniquely defined by (j5.4.8[) . 

Substituting (15.4.80 into (|5.4.5|) gives 

sn2 (y) = t^ + ^{s^ + l) - ^ ^(2*2 + ^2 -1)2 + 4^2. (5.4.9) 

Here we need only notice that dn(x) is strictly positive to see that sgn(sn(y)) = sgn(t) 
giving 

sn{y) = sgn(t)^t2 + 1 (^2 _^ 1) _ iY^(2i2 + - if + 4s2. (5.4.10) 

The fact that sn is bijective on [—yo,yo] completes the proof. □ 
5.4.2 G and H are semi-conjugate 

The cornerstone of our proof that Q : A ^ A is Bernoulli is the result that it is 
semi-conjugate to H. The formal statement follows and will be proven in this section. 

Proposition 5.4.7. The identity E o H = Q o E holds on R. 

We first show that E o F = ^ o E for (x, y) G P. We deal separately with the cases 
(x,y) G Pq and (x,y) G Pi, the former being trivial: 

Lemma 5.4.8. The identity E o F = ^ o E holds on Pq. 

Proof. By definition $ : A+ is given hy E o T o E~^, where by E~^ we mean 

the inverse to the restriction of to Pq- T : C ^ C is the linear twist map defined 
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in Section ri.2.21 so it suffices to show that given {x,y) G x [—yo,yo] then T{x,y) = 
F{x,y). This is immediate from the respective definitions. □ 



Proving that EoF = ^oE on Pi will require a little more work. Let J : — > be 
defined by J(x, y) = {2K — x, 2K — y). It is easy to establish that J is a diffeomorphism, 
is its own inverse and that J{Pq) = -Pi- S 

Lemma 5.4.9. The identity E o J = E holds on T^. 

Proof. The proof is an elementary application of Proposition 15.4.1] E{J{x,y)) is given 

by 

{sn{2K - 2;)dn(2i^ - y),cn{2K - x)cn{2K - y),dn{2K - x)sn{2K - y)) 
= (sn(x)dn(y), — cn(x) • — cn(y), dn(x)sn(y)) 
= E{x,y). 

□ 

Lemma 5.4.10. The identity J o F = F o J holds on Pq. 

Proof. The proof is straight-forward calculation. Let {x,y) £ Pq, so J{x,y) £ Pi, then: 

FoJ{x,y) = F{2K -x,2K -y) 

= {2K-x + {AKyo - 4Ky)/2yo, 2K - y) 
= (2K-x + 2K- 2Ky/yQ, 2K - y) 

and conversely 

JoF{x,y) = J{x + AK{y + yo)/2yo,y) 

= {2K-x- 2Ky/yQ - 2K, 2K - y) . 

J is sometimes referred to as an involution and reflects points through the origin in We 
are grateful to Prof. Robert MacKay for pointing out to us the fact that T'^ with all pairs of points 
{(x, y), J(a;, J/)} identified is topologically equivalent to S^. Our map E provides an explicit means by 
which one might relate a topology on to the quotient topology on . 
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The two are equal because 2K = —2K in S^. □ 

It is now simple to prove the following: 
Lemma 5.4.11. The identity E o F = ^ o E holds on Pi. 

Proof. Let {x,y) £ Pi, let {x,y) = J{x,y) and observe that {x,y) £ Pq. We must show 
that E o F{x, y) = ^ o E{x, y). Using Lemma 15.4.91 we have 

^oE{x,y) = ^oE{x,y), (5.4.11) 

and using Lemma 15.4.101 followed by Lemma 15.4.91 we have 

E o F{x, y) = EoJo F{x, y) = E o F{x, y). (5.4.12) 

Lemma 15.4.81 says that the expressions (|5.4.1ip and (|5.4.12p are equal. □ 

In very similar fashion it can be shown that E o G = T o E holds on Q. The 
semi-conjugacy follows immediately: 



Proof of Proposition \5.4-T\ Let {x,y) G i?, then E o G o F{x,y) = T o E o F{x,y) 



T o^o E{x,y). □ 
5.4.3 Proof of Theorem [TTS^l 



Ornstein (197ll ) and essentially follows 



Our main result is a consequence of a theorem off* 
from the semi-conjugacy just established, although we will need some supplementary 
results. We need to discuss measure- theoretic factors and we begin by reviewing a few 



definitions; we have taken these from 



Katok and Hasselblatt (1993 ) 



A measure space (X, A4,z^) with finite measure v is called a Lebesgue space if it is 
isomorphic to the union of [0, 1] with Lebesgue measure, with at most countably many 
points of positive measure. 

124 



For measure-preserving transformations Ti : Xi — > Xi and T2 : X2 X2 of 
Lebesgue spaces (XijUi) and (X2,f2) respectively, we say that T2 is a metric factor 
of Ti, or from now on just a factor of Ti, if there exists a measure preserving map 
9 : Xi ^ X2 such that 

BifVi = V2 and T2 o 6 = 8 o Ti. 

The notation ^^i^i denotes the pushforward measure on X2 obtained from the measure 
on Xi: for each measurable set B C X2 this is defined by 



We will use the following: 



Theorem 5.4.12 (jOrnstein (197ll )). A factor of a Bernoulli map is itself Bernoulli. 



We prove some results concerning Q and {A,fl), where ft = E^,^. Let M. denote the 
(T-algebra of Lebesgue-measurable subsets of R and let 



M = {B(ZA: E-\B) G M}, 



where E ^(B ) denotes the pre-image of i? C ^ with respect to E. It is known (see 



Rudin (19871 ). Theorem 1.12, p. 13) that Ai is a fi-algebra of subsets of A. 

Recall that /i denotes the Lebesgue measure on R. Let jl : M. ^ [0,1] be given by 

fL = E^,fl. 

Proposition 5.4.13. {A,fl) is a Lebesgue space. 

Proof. The function /i, defined on u-algebra M, takes its range in [0, 1]. We show that 
it is countably additive and thus a measure. Let {Bi : z G N} be a disjoint, countable 
collection of members of Ai, so {E^^(Bi) : i G N} is a disjoint, countable collection 
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of members of A4. By the countable additivity of fi we have fi [^U£i ^i-^i) 
X^i^i l^i^^^i^i))- The result then follows from the observation: 

\jr=iE~'m = U=i {{u,v)eR: E{u,v) e B,] 
= {{u,v) e R : E{u,v) e UZiBi] 



Next, we appeal to the result that any Borel probabilit y measure on a separable, lo- 



Katok and Hasselblatt (19951 ). 



cally compact H aus dor ff space defines a Lebesgue space (see 
Theorem A. 6. 7, p. 734). It is not difficult to see that (A, fl) satisfies these criteria: A 
is a separable Hausdorff space because has these properties and by the Heine-Borel 
theorem it is compact and hence locally compact. It is obvious that fi{A) = 1; to see 
that Ai contains all Borel subsets of A, let C A be open, then E^^{B) C i? is open 
(because E is continuous) and so E~^{B) € M^. By definition B e M.. □ 



Katok and Hasselblatt (19951 ) shows equally that (i?, /i) 



The result we quoted from 
is a Lebesgue space. 

Lemma 5.4.14. Q : A ^ A preserves fl, i.e. if B £ Ai then 

m{B))=KB). (5.4.13) 

Proof. The proof involves a little manipulation of identities we have established. Let 
B eMand let B = E-^{B) e M. It follows that E{B) = B and 

B = E-'^ o E{B). (5.4.14) 

From Lemma [5.4.7l it follows that H^^oE^^ = E^^oQ^^ (note that these are pre-images 
rather than functions). Using this fact and (|5.4.14p we deduce that 
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B = E~^ o E{B) = E-^ o o e o E{B) = H'^ o E'^ o E o H{B) 



I.e. 



H{B) = E-^ oEoH{B). (5.4.15) 
By definition we have 

fi{B) = iJi[E-\B))=iJi{B) (5.4.16) 
and using (j5.4.14p followed by (j5.4.15p we have 

fl (e{B)^ =fioE~^oQo E{B) = fioE~^oEo H{B) = n {H{B)) . (5.4.17) 

The expressions ()5.4.16p and (|5.4.17p are equal because H preserves /Lt. □ 

Our main result follows easily. 

Proof of theorem 1^.3.31 Propositions 15.4.71 and 15.4.131 and Lemma 15.4.141 show that 

is a metric factor of H. Theorem 15.4.121 completes the proof. □ 



1971 



paper in which Theorem 15.4. 121 is 



We end with a remark regarding Ornstein's 
established. Ornstein defines a factor of a Bernoulli shift to be the restriction thereof to 
an invariant sub-ci-algebra; we demonstrate briefly that the (more common) definition 
we have taken is equivalent. Indeed, if we let M' C M consist of precisely those elements 
B £ Ai for which B = E^^ o E{B) and let /u' be the restriction of /u to M' then it 
follows from the results of this section that E : (R, fi') — > {A, ji) is an isomorphism and 
thus the result. 
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6 Summary and outlook 



We finish by surveying the results we have estabhshed and discussing some strengths and 
weaknesses of our methods. We consider the directions in which productive future work 
might be undertaken, either as a direct consequence of the present work or otherwise. 
At the end of the chapter we will propose two conjectures which, we believe, would be 
an excellent starting point for anyone who wished to generalise our methods to the class 
of abstract linked-twist maps we have defined. 

6.1 Summary 

We make some comments about the results we have obtained. 

6.1.1 A topological horseshoe in the toral linked-twist map 

In Chapter [3] we established the existence of a topological horseshoe in th e toral linked- 



Devanev fl978l ) 



twist map defined in Section [1.3. 11 Our method was inspired by the work of[ 
who constructed such a horseshoe in the planar linked-twist map defined in Section ri.3.21 
We observed an interesting difference between the two constructions. The conjugacy 
constructed by Devaney is with a sub-shift of finite type, conversely ours is with full shift 
on N symbols. 

We believe that this can be explained by the fact that the planar linked-twist map 
has two distinct intersection regions, so that in this case the invariant Cantor set A is 
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split between the two. We conjecture, based on our result, that there is an invariant 
Cantor set A+ C S+ for the planar map on which the dynamics, as in our example, are 
conjugate to a full shift on the space of symbol sequences. This set may be constructed 
by considering only those points that land in on each iteration, and excluding the 
other points from Devaney's construction, so in fact it is a proper subset of his invariant 
set. 

Devaney's invariant set has a richer structure than the one that we have constructed, 
and we believe this to be a consequence of it containing as proper subsets a wealth of 
other invariant sets whereby there is some restriction on which of S-t points return at 
any given time. An investigation of this structure would certainly be an interesting 
exercise in its own right, although it is not clear that it would yield any conclusions 
that might help us to better understand the dynamics on a set of full measure. 

6.1.2 The Bernoulli planar linked-twist 

The shortcomings of our result are clear: we have established the Bernoulli property 
for a planar linked-twist map composed of the embeddings of linear twists, but we have 
had to be explicit about the sizes of the annuli, taking the inner annuli to have size 
ro = 2 and the outer annuli to have size ri = \fl . It is clear where these restrictions 
were required so let us look at this a little more closely. 

Key to our proof was Proposition 14.4.1] which states that, with annuli of the sizes 
specified, then DF preserves the tangent cone C = {(u, v) : uv > 0} and DF^^ preserves 
the tangent cone C = {{u,v) : uv < 0}. We do not intend to repeat all of the details here, 
but an example illustrates the difficulty. We were required to show (see equation ()4.4.3p ') 
that y) > , or to give the full expression in terms of the function ip, that 



Diip{x,tp ^{x,y) + c{x - ro)) 

+D2tp {x,ij~^{x,y) + c{x - ro)) Di'4)~^{x,y) + 



2lT 



>0, 
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for each pair {x,y) S [^0,7-1] x [0,7r]. 

We took a rather crude approach to this and sought to bound each of the terms 
Diip{-,-), D2^{-,-) and Di7p~^{-,-) individually. We observed a lower bound of for 
Diip which (from the proof of Lemma [4. 4. 3p seems optimal, but it is quite possible that 
none of the other bounds established are optimal. Of course, even if we were to obtain 
optimal bounds on each of the three derivatives individually this would not necessarily 
give us optimal bounds for Dif^{x,y). 

Given our crude approach to this problem it is perhaps remarkable that it works for 
any system at all. The fact that we were able to find (and, we should add, with relative 
ease) choices of rg and ri for which the problem is tractable could be interpreted as 
evidence that the inequalities in fact hold for a much wider choice of annulus size. 

If one is motivated to use our method to prove the Bernoulli property for some range 
of ro and ri values then a more sophisticated approach to these inequalities is imperative. 
Plotting Dif^{x, y) over the required domain would be a good start, although even this 
is non-trivial as one needs a package with sufficiently good programming and graphical 
capabilities, due to the nature of the functions il) and Such a plot might suggest a 
way to partition the domain so that tighter bounds can be established on each partition 
element; this would seem to entail a great deal of work however. 

We conclude by discussing how far one might hope to develop the method we have 
introduced. The ultimate ambition would be to give a 'complete description' of the 
possible dynamics and this, perhaps, might consist of a large open sets of parameter 
values for ro,ri where the Bernoulli property is established, and a complementary set 
on which it is shown not to occur. 



Sturman et al. (20061 ) provide a number of plots showing numerical simulations of 
(co-twisting) planar linked-twist maps, some of which appear to exhibit good mixing 
properties and others which do not. Recall that one of our initial assumptions was 
of transversality, which we believe to be related to our ability to construct the new 
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coordinates (we will say more on this shortly). It would seem from the simulations that 
transversality is not a prerequisite for good mixing, and hence the method we have 
proposed cannot be expected to provide such a complete description as we have asked 
for. Of course, this leaves open the possibility that transversality is sufficient for a 
linked-twist map to be Bernoulli. 

6.1.3 The Bernoulli linked-twist map on the sphere 

It is certainly interesting that we have been able to construct so directly a semi- 
conjugacy between a map on the torus and a map on the sphere. This is an immediate 
consequence of the coordinate system we have used. The coordinate transformation 
would perhaps be of interest to the wider mathematical community given its relatively 
clean expression and the orthogonal coordinate system it provides. The most interest- 
ing development from a dynamical systems perspective might be to use it to construct 
further examples of Be rnoul li maps on full measure subsets of the sphere, by a method 



analogous to Katok's (|l979l ). Recall that the starting point for his construction is a 



hyperbolic toral automorphism with certain points fixed. 

Perhaps the strength of the method is also its weakness; it is quite specialised and so 
it is difficult to see how one might hope to generalise it in order to obtain other results, or 
indeed what those other results migh t be. Neverthele ss it afforded us the opportunity to 
use some techniques (the theorem oi[ 



Ornstein (197ll )) that perhaps otherwise we would 



not have discovered. 



6.2 Ideas for further work 

We conclude by looking at two ways in which one might build upon the results we have 
established. 
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6.2.1 Decay of correlations 

We have mentioned many applications for which certain hnked-twist maps provide a 
natural model and thereby a means to understand or to predict the behaviour to some 
extent. This is of particular importance within the nanoscale devices we have mentioned 
such as the DNA microarray discussed in Section [2.4. 11 because the alternative trial and 
error approach to their design is prohibitively costly. We have mentioned that there is a 
degree of convergence between those questions that are interesting from a mathematical 
perspective and those that are interesting from an applications perspective. 

The strength of mixing is one such question. We have established for two different 
types of linked-twist map that the Bernoulli property is satisfied. The implication for 
systems whose dynamics are well approximated (in some sense) by these maps is that 
they should be expected to mix initial conditions thoroughly. 

The rate of mixing is another such question which we shall discuss briefiy now. 
Consider for a moment the cornerstone of our proofs; we have spelled this out previously 
but we re-iterate it now. There is a region of positive area (which we have labeled S 
in each case) to which almost every point returns an infinite number of times. The 
hyperbolicity, which is responsible for the separation of nearby trajectories and thus 
the strong mixing, is inextricably linked to this behaviour. 

Now consider the size (i.e. the measure) of this region relative to the size of the 
whole manifold A. Bernoulli systems automatically satisfy the strong mixing property 
(defined in Section I2.1.ip which says that for a measurable set B G A having positive 
measure 

y M(e"(i?)ns) 

where we have used the invertability of 0. We might interpret this as saying that the 
asymptotic proportion of i? in S is proportional to the size of S. Given the relationship 
between returns to S and separation of nearby trajectories, one might conjecture that 
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the greater the size of S, the more 'chaotic' the system is in some sense and, importantly 
for apphcations, the faster the phase space becomes mixed. 

The concept of the decay of correlations is the correct mathematical formalism 
within which to phrase such questions. In essence the idea is to look at the rate of 
converge nce in the ergo dic theorem. We don't provide a formal definition but direct the 



reader to 



Baladi (200ll ) for further details. 



The seminal work in r ecent years on . the decay of correlations in dynamical systems 



with some hyperbolicity is lYoung (19981 ) . She establishes a framework for studying this 
decay in a class of systems she characterises as having 'regular returns to sets with good 
hyperbolic properties '. In this context we see the importance of the perspective we have 
taken in analysing linked-twist maps and why we are hopeful that this approach will 
prove useful in future endeavours. 

6.2.2 Ergodic properties of abstract linked-twist maps 

Our definion of an abstract linked-twist map invites the question of its ergodic prop- 
erties. We discuss briefly how the proofs of general results along these lines might be 
attempted using the techniques introduced in this work. We stress that this should 
not be considered a 'work in progress'; rather these are merely preliminary comments 
which we hope may be of inspiration to anyone inclined to persue results along these 
lines, and may a t least be of some interest to other readers. It is worth remarking that 



Przvtvcki (19811 ) appears to have attempted results along these lines. To what extent 
he has acheived these ambitions is not entirely clear to us. 

Recall our definition of an abstract linked-twist map on a smooth manifold M of 
dimension 2: for i = 1,2 let Cj be a cylinder, let Ei : Ci ^ A C M he an embedding, let 
Ti : Cj ^ Cj be a twist map on Cj and let E'j o Tj o E^^ be a twist map on Ei{Ci) C R. 
If Ei{Ci) and i?2(C'2) are transversal then the composition 
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Q = E20T20 E^^ oEioTioE- 



.-1 



(6.2.1) 



is called a linked- twist map. 

In analysing the linked-twist map in the plane it was crucial that intersection regions 
S-t (that is, the connected components of £'i(Ci) H £'2 ((^2)) could each be expressed in 
new coordinates in which it was the Cartesian product of two intervals (a 'square', in 
fact). (In the case of the sphere this followed immediately from our choice of coordi- 
nates.) We conjecture the following: 

Conjecture. Transversality of the embedded cones Ei{Ci) and E2{C2) is a sufficient 
condition for the existence of local coordinates in which each connected component of 
-Ei(Ci) n £'2(^2) is a Cartesian product of two intervals. 



Figure 6.1: An example of an abstract linked-twist map. Part (a) shows the view parallel 
to the y-axis and part (b) parallel to the x-axis. The manifold A (shaded) consists of 
two cylinders Ci and C2 embedded into with embeddings Ei and E2 respectively. 

We illustrate an example of the situation we have in mind in Figure 16.11 Q Part (a) 
shows a projection onto the xz-plane, whereas part (b) shows a projection onto the 

^We thank Prof. Jens Marklof for suggesting this map to us. It is a linked-twist map defined on S'^ 
and, as opposed to the map studied in Chapter [S] points moving under a twist map Ei oTi o Ei do so 
in a plane of constant x 01 z coordinate (where {x, y, z) are Cartesians in M^). Moreso than the map of 
Chapter [5] this resembles the motion of a 'top' undergoing precession. It is therefore possible that this 
map might embody the essence of certain quantum chaotic motion. 




(a) 



-Xq Xq 



(b) 



-2/0 2/0 
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yz-plane. The two embedded cylinders Ei{Ci) and £'2(^2) are bounded by lines of 
constant z and x coordinate respectively and their union is denoted by ^ C S'^. 

In Figure 16.21 we show the cylinders themselves with the pre- images of the intersec- 
tion regions shaded. Let (si, ii) G x /i give coordinates on Ci and let (52,^2) £ x I2 
give coordinates on C2. 

§1 

s^^ 

Figure 6.2: The two cylinders with pre-images of the intersection regions shaded. 

Let ^ denote Lebesgue measure on A; it follows by the usual arguments (assuming 
that the Tj are sufficiently well behaved) that fi-a.e. point returns infinitely many times 
to Define the first return map Qy. '■ ^^^{T,) E^^{Yi) by 

= 0E20 o E^^ oEio 

where n and m are positive integers satisfying the usual criteria. Finally, define the 
usual tangent cone C consisting of the open first and third quadrants. If D {E2^ o -Ei) 
and D [E^^ o E2) preserve and expand the cone C then DQ^, will also. In this case we 
conjecture that Q has the Bernoulli property. 

If this isn't the case (a situation analogous to the planar linked-twist, where a 
larger cone U was preserved but C was not) then we might still be able to proceed as 
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before. Using the ideas of Chapter H] we can construct new coordinates (if the previous 
conjecture holds) on A, which are equal to (11,12) for a point {u,v) £ T,± such that 
Ej^{u,v) = {ij,Sj), j = 1,2. The coordinate transformations we have mentioned can 
be expressed in terms of the embeddings i?i , £'2 ■ 

Conjecture. We can establish sufficient criteria for Q to have the Bernoulli property. 
These criteria consist of a pair of inequalities involving only the derivatives of E^^ and 
Ti,i = 1,2. 
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